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ABSTRACT 

Context. Extragalactic observations allow the study of molecular chemistry and excitation under physical conditions which may differ 
greatly from what found in the Milky Way. The compact, obscured nuclei (CON) of luminous infrared galaxies (LIRG) combine large 
molecular columns with intense infrared (IR), ultra-violet (UV) and X- radiation and represent ideal laboratories to study the chemistry 
of the interstellar medium (ISM) under extreme conditions. 

Aims. To obtain for the first time a multi-band spectral scan of a LIRG, in order to derive the molecular abundances and excitation, to 
be compared to other Galactic and extragalactic environments. 

Methods. We obtained an ALMA Cycle 0 spectral scan of the dusty LIRG NGC 4418, spanning a total of 70.7 GHz in bands 3, 6 , 
and 7. We use a combined local thermal equilibrium (LTE) and non-LTE (NLTE) fit of the spectrum in order to identify the molecular 
species and derive column densities and excitation temperatures. We derive molecular abundances and compare them with other 
Galactic and extragalactic sources by means of a principal component analysis. 

Results. We detect 317 emission lines from a total of 45 molecular species, including 15 isotopic substitutions and six vibrationally 
excited variants. Our LTE/NLTE fit find kinetic temperatures from 20 to 350 K, and densities between 10^ and 10^ cm“^. The spectrum 
is dominated by vibrationally excited HC 3 N, HCN, and HNC, with vibrational temperatures from 300 to 450 K. We find that the 
chemistry of NCG 4418 is characterized by high abundances of HC 3 N, SiO, H 2 S, and c-HCCCH and a low CH 3 OH abundance. A 
principal component analysis shows that NGC 4418 and Arp 220 share very similar molecular abundances and excitation, which 
clearly set them apart from other Galactic and extragalactic environments. 

Conclusions. Our spectral scan confirms that the chemical complexity in the nucleus of NGC 4418 is one of the highest ever observed 
outside our Galaxy. The similar molecular abundances observed towards NCG 4418 and Arp 220 are consistent with a hot gas-phase 
chemistry, with the relative abundances of SiO and CH 3 OH being regulated by shocks and X-ray driven dissociation. The bright 
emission from vibrationally excited species confirms the presence of a compact IR source, with an effective diameter <5 pc and 
brightness temperatures >350 K. The molecular abundances and the vibrationally excited spectrum are consistent with a young 
AGN/starburst system. We suggest that NGC 4418 may be a template for a new kind of chemistry and excitation, typical of compact 
obscured nuclei (CON). Because of the narrow line widths and bright molecular emission, NGC 4418 is the ideal target for further 
studies of the chemistry in CONs. 
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1. Introduction 

Extragalactic chemistry is a field that is quickly expanding lead¬ 
ing to new, powerful diagnostic tools for the star-for ming and ac¬ 
tive galactic nuclei (AGN) activity in galax ies (e.g., iMeier et al.l 
l2014t IViti et all I2014I : iMartm et al.l l2015h . The extreme envi¬ 
ronments found in some extragalactic objects provide the op¬ 
portunity of studying the properties of the interstellar medium 
(ISM) beyond the typical conditions found in the Milky way. 
Shocks, stellar- and AGN radiation, dust shielding, and cosmic 


rays strongly impact the chemistry and excitation of the molec¬ 
ular ISM. Establishing the chemical and physical conditions of 
the molecular gas becomes a particularly important identifica¬ 
tion tool when the activity itself is buried in dust. 

A clear case is offered by the compact obscured nuclei 
(CON) of IR-lu minous (LIRGs) and ult raluminous galaxies 
(ULIRGS, e.g., ISanders & Mirabel 1 1996h . These galaxies ra¬ 
diate most of their energy as thermal dust emission in the in¬ 
frared and constitute the dominant population among the most 
luminous extragalactic objects. Observations at mid-IR and mil- 
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limeter wavelengths suggest that they may play a crucial role in 
galaxy evolution, represe nting the early obscured stages of star- 
burst galaxies a nd AG Ns (ISnoon et ahl 20071: Aalto et ^l2007al : 
ISakamoto et al.l 1201ICostagliola et al. 2013h . Because of the 
high amounts of gas and dust (N(H 2 )> 10 ^'* cm“^) and the ex¬ 
treme obscuration (Av>100 mag) the nature of the activity can¬ 
not be reveled by standard optical, IR or even X-ray observa¬ 
tions. This makes it essential to identify new diagnostic tools in 
the mm/submm window, which can probe deeper in the column 
of dust. 


Until recently, research in this field has followed mainly 
two paths: molecular line-ratio surveys, and single-band spec¬ 
tral scans. Line-ratios surveys have mostly focused on the 
brightest transitions of dense-gas tracers such as HCN, HNC, 
and HCO^, trying to compare the observed values to chemi- 


inated regions (XDRs), and hot-cores (e.g 

Imanishi et al.l 

20071 Loenen et al. 12007: Krios et alJ 

20081 

iBaan et al.l 12008: 

Costagliola et al.ll201 it Viti et al.ll2014t) 

However, the line ratios 


of such bright tracers show only small variations, of the order 
of a few, and their theoretical interpretation is still controver¬ 
sial (e.g..lKohno et al. 200 lH Aalto et al.ll2007bl:lMeiierink et^ 


1200'^ ICostagliola et alj201 Ih . Since these molecules are ubiqui¬ 

tous in the dense ISM and are excited under a variety of physical 
conditions, the observed trends may be due to opacity and exci¬ 
tation effects, rather than a different chemistry. 

Unbiased spectral scans are powerful tools to find new trac¬ 
ers, more sensitive to the physical conditions of the gas. How¬ 
ever, wide-hand molecular line surveys have been carried out 
only for a very li mited number of extragalactic objects (e.g ., 

iMartm et al.ll2^l201 UfMuller et al.ll201 iblAladro et^l20iai . 


While these obseryations haye made a leap forward in the study 
of the chemical complexity of galaxies, most of them haye been 
limited by the narrow frequency range coyered, which does not 
allow a complete analysis of the molecular excitation. In order 
to compare chemical models with the obseryations, a better han¬ 
dle on the abundances of the molecular species is needed. This 
can be achieyed only by obserying multiple transitions in differ¬ 
ent frequency bands, sampling a wide range of excitation states. 
Here we present the first multi-band spectral scan of the ob¬ 
scured luminous infrared galaxy NGC 4418 obtained with the 
Atacama Large Millimeter/submillimeter Array (ALMA) obser- 
yatorjQ. 


1.1. NGC 4418 : The prototypical obscured LIRG 

The LIRG NGC 4418 (Lir = 10" Lq) has the optical morphol¬ 
ogy of an early-type spiral and was first detected as a bright 
source by the IRAS satellite. Lying at D-34 Mpc, the galaxy is 
part of an interacting pair, with the companion being an irregular 
blue galaxy roughly two arc-minutes (16 kpc) to the southeast. 

The optical spe ctrum of NGC 4418 has been described by 
iRoche et~51 (Il986h as unremarkable, with only faint emission 
from S[II] and Hq, on a strong continuum, and it does not hint 
at the presence o f a bright IR source . This is explained by mid- 
IR obseryations (ISnoon et al.ll2001at iRoche et'^ll986h . which 
show a deep silicate absorption at 10 pm, one of the deepest eyer 
obseryed, corresponding to roughly 100 magnitudes of yisual ex¬ 
tinction. Whateyer is powering the strong IR flux of NGC 4418, 
it must be hidden behind a thick layer of dust, which makes it 
undetectable at optical wayelengths. 


* http://www.almaobservatory.org/ 


The high IR luminosity requires a compact starburst or an 
AGN to heat up a large column of dust. However, because of 
the high obscuration, the direct investigation of the nuclear re¬ 
gion is a challenging task, and the nature of the energy source 
is still unclear. Observ ations with the Chandra X-ray satellite by 
iMaiolino et al.l (l2003l) show a flat hard X-ray spectrum, which 
would be an indication of a Compton-thick AGN, but the photon 
statistics are too limited to be conclusive. The absence of a clear 
X-ray signature may imply either that the galaxy is powered by 
star formation alone, or that the obscuring column is so high that 
most of the X-ray emission cannot escape its nuclear region. 

Additional evidence of nuclear activity in NGC 4418 comes 
from near- and mid-IR observations. High-r esolution HST near- 
infrared and Keck mid-infrared images by lEvans et aP (l2003h 
show that the galaxy has a 100-200 pc optically thick core, with 
a high IR surface brightness. The observed spectral energy dis¬ 
tribution implies a dust temperature of 85 K, which, when com¬ 
pared with the total IR flux, implies the presence of an optically 

thick source of no more than 70 pc across. _ 

_ ^ two recent papers (ISakamoto et al.l2013l : ICostagliola et"^ 

12013h we used the Submillimeter Array (SMA) and the Multi- 
Element Radio Linked Interferometer Network (MERLIN) to 
directly probe the nucleus of NGC 4418 at mm and radio wave¬ 
lengths with sub-arcsecond resolution. These studies confirm the 
existence of a ~20 pc ( 0 ." 1 ) hot dusty core, with a bolometric lu¬ 
minosity of about 10*' Lq, which accounts for most of the galaxy 
luminosity. The high luminosity-to-mass ratio {L/M ^^500 L© 
Mg') and luminosity surface density Lg pc“^) are con¬ 

sistent with a Compton-thick AGN to be the main luminosity 
source. Alternatively, an extreme (SER^^lOO Mg yr“*), young 
(<5 Myr), compact starburst could also reproduce the properties 
of the inner core. 

Combined observations with MERLIN and the Euro- 
pean Very Lo n g Bas eline Interferometry Network (EVN) by 
IVarenius et al.l (l2014l) show that 50% of the total radio emis¬ 
sion at 5 GHz originates from eight compact sources in the inner 
0."2 of the nucleus. The four brightest sources lie inside a 0."1 
circle centered on the position of the 860 pm continuum from 
ISakamoto et al.l (1201.31) . These sources have an average surface 
brightness which is close to the limit of what can be produced 
by well-mixed thermal/n on-thermal emi s sion f rom any surface 
density of star formation. IVarenius et al.l (12014l) suggest that the 
radio emission could be explained by super star clusters with in¬ 
tense star formation with some contribu tion from an AGN. 

He rschel PACS observations by iGonzalez-Alfonso et al.l 
( 2012h. and our combined MERLIN/SMA study 
( Costagliola et al.l 12013h reveal the presence of redshifted 
OH and HI absorption, which is interpreted as the signature 
of a molecular inflow. The presence of a molecular outflow 
is also suggested by a U-shaped red optical f eature along the 
north western semi-minor axis of the galaxy (ISakamoto et al.l 
l2013l) . 


NG C 4418 was f irst sho wn to have a rich molecular chem¬ 
istry by jAulto_etayj2007^ and the high abundance of HC 3 N 
(>10“^, ICostagliola & Aaltal2010h is not reproduced by mod- 
els of X-ray-dominate d chemistry expected in an AGN scenario 
(iMeiierink et aLll^OOTll . Together with a low HCO^/HCN 7=1-0 
line ratio, bright HC 3 N is instead reminiscent of line emission 
toward Galactic hot-cores, i.e., regions of dense, warm, shielded 
gas around young stars. This has led some authors to propose 
that NGC 4418 may be a very young starburst, where the star- 
forming regions are still embedded in large amounts of dus t 

(iRoussel et al.ll200^ I Aalto et ^l2007at ICostagliola et al.ir201 111 . 

This scenario of a nascent starburst would be consistent with the 
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galaxy being radio-deficient. How ever, recent results in chemical 
modeling bv iHarada et alJ (1201 3h suggest that substantial HC 3 N 
abundances may be maintained near an AGN. 

The combination of a bright, compact IR source and a 
high molecular column produce an incredibly rich vibrationally 
excited spectrum. Vibrationally e xcited HCN and HNC were 
first detected in this s ource by ISakamoto et alJ (l 2010 h and 
ICostagliola et al.l (1201 3h . and the galaxy shows very bright 
emission from yibrationally e xcited HC 3 N (ICostagliola & Aalt^ 
l 2010 HCostagliola et al.ll2013h . 

The yelocity dispersion of the molecular emission lines 
emerging from the core of NGC 4418 is small, of the order 
of ~100 kms“^ This makes the identification of spectral fea- 
tures much more re liable than in other (U)LIRGS (e.g., Arp 220, 
iMartm et al.ll201 1|) . The combination of bright molecular emis¬ 
sion and narrow lines makes it the ideal target to study molecular 
chemistry and excitation in compact obscured nuclei. 

2. Observations 

The obseryations were undertaken during the ALMA Early Sci¬ 
ence Cycle 0 phase, between April 2012 and January 2013. The 
target NGC 4418 was obseryed in 11 tunings out of 38 proposed 
(see Table [T]i. Each tuning consists in short integrations of ~5 
minutes on source, plus obseryations of the bright radio quasars 
3C 273 or 3C 279 and Titan for absolute flux calibration (except 
for 3 dates, for which we used 3C 273). Eor each tuning, four dif¬ 
ferent 1.875 GHz-wide spectral windows were set, each counting 
3840 indiyidual channels separated by 0.488 MHz. 

In order to account for possible edge instabilities and offsets, 
contiguous spectral windows were oyerlappedby 250 MHz. As a 
result, each tuning coyers a total span of 7 GHz in sky frequency, 
including upper and lower sideband. The obseryed tunings coyer 
28.7 GHz in band 3, 28 GHz in band 6 , and 14 GHz in band 7, 
corresponding to the 90%, 43%, and 14% of each band, and a 
total coyerage of 70.7 GHz. 

The Cycle 0 array configuration was composed of 19 to 25 
antennas, with baseline separations from 20 to 480 meters. The 
resulting angular resolution is of ~ 2 ", ~ 1 ", and ~ 0".8 for bands 

3, 6 , and 7 respectiyely. 

The data calibration was done within the CASA0 package 
following a standard procedure. We expect an absolute flux ac¬ 
curacy of the order of ~10% in band 3 and ~20% in bands 6 and 
7. 

2. 1. UV-fit of the spectrum 

Pre yious interferomet r ic obse ryations by lSakamoto et al.) (l2013l) 
and ICostagliola et al.l (l2013l) show that 70% of the molecular 
emission in NGC 4418 emerges from a region of less than 0".4 in 
diameter, we therefore expected all the emission to be included 
in our synthesized beam. We yerified this assumption by fitting 
the source size of the emission of the HNC 1-0 and HNC 3-2 
lines in band 3 and 6 . We performed both a fit of the yisibilities 
and a 2D Gaussian fit of the integrated emission in the image 
domain. We found the emission to be unresolyed, with upper 
limits to the source size of 0".4 and 0".2 for band 3 and band 6 , 
respectiyely. 

Giyen the point-like nature of the emission, we chose not to 
image and clean the whole dataset, but to extract the spectrum 
from the yisibilities. We first obtained the position of the contin¬ 
uum peak by ayeraging line-free channels in each tuning. Then 

^ http://casa.nrao.edu/ 


we fitted the interferometric spectral yisibilities of each spectral 
window assuming a point-source model, with a fixed position 
but free amplitude. Th is last step was performe d using the CASA 
uvmultifit routine (iMartf-Vidal et al.ll2014l) . 

The extracted channels were then interpolated to a common 
resolution of 20 km s ' across the whole frequency range, and all 
spectral windows were merged to create a single spectrum. Small 
offsets between contiguous spectral windows were eliminated by 
ayeraging the channels in the 250 MHz oyerlap region. These 
offsets were not larger than 10 % in either band, confirming the 
good quality of the data calibration. The measured rms of the 
spectrum is 2 mJy in band 3, 6 mJy in band 6 , and 4 mJy in band 
7. The total extracted spectrum is shown in Eig.[T] 

In Eig. lA.lOl we compare the extracted ALMA spectrum with 
IRAM 30 m single-dish obseryations with beam sizes of 27" and 
9" at 3 and 1 mm, respectiyely. The flux densities measured for 
the two datasets are within the errors, which is consistent with all 
the molecular emission to be contained inside the ALMA beam. 

A comparison with the HNC 3-2 and HC 3 N 30-29 flux 
densities detected with the SMA in the extended configuration 
(ICostagliola et al.ll2013[) is reported in Eig. lA.llI The integrated 
flux density of the different obseryations is reported in Table lA. 8 l 
We find that between 70% and 90% of the total flux is contained 
by the 0".4 synthesized beam of the SMA. 

3. Molecular emission 

The mm-waye spectrum of NGC 4418 shows rich molecular 
emission in all three ALMA bands. Despite the narrow line 
width of ~120 km s ', line confusion is reached after short in¬ 
tegrations of ~5 minutes in bands 6 and 7. 

3.1. Removing the continuum 

Remoying the continuum from the spectrum in band 6 and 7 is 
complicated by line crowding. In order to fit a continuum leyel, 
we assumed the regions in bands 6 and 7 with the lowest emis¬ 
sion in the band to be line-free. We then ayeraged the channels 
to obtain band-ayeraged yalues for the continuum. The deriyed 
yalues are 9+1 mJy, 56+6 mJy, and 85+9 mJy at 98, 244, and 
284 GHz, respectiyely. 

In Eig. 13 we compare the der iyed continuum yalues with 
data from SMA obse ryations by ISakamoto et al.l (120101) and 
ISakamoto et al.l (l2013h . Allowing for a calibration uncertainty 
of 10 %, the two datasets are consistent. 

We find that the emission is best fitted by two power laws 
(Sy cc v“), one between 97 and 300 GHz with a = 2+0.4, and 
one at frequencies higher than 270 GHz with a - 3+0.6. If we 
assume a uniform slab geometry, the dust spectral index in the 
Rayleigh-Jeans regime can be written as a -2+/3T{e^ - 1)“*, 
where /J is the the index of power-law frequency dependence 
of the dust opacity and t the opacity of the emission at the 
wayelength where a is measured. Eor a typical yS=L5-2 (e.g., 
iLisenfeld et al.ll2000l) . an a -3 requires opacities of the order 
of un ity. This result is in ap eement with preyious obseryations 
(e.g., ISakamoto et al.l2013l) which required an optically thick IR 
source to fit the galaxy’s spectral energy distribution. The yalue 
of O' = 2+0.4 at lower frequencies is difficult to explain as an 
opacity effect because it would require r > 10 and one would 
expect opacity to increase with frequency. A better explanation 
may be that free-free or cold dust emission is contributing sig¬ 
nificantly to the continuum in band 3. 

The best power-law fit to the ALMA continuum points oyer- 
estimates the emission at 1 mm by more than 20 %, too much for 
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Table 1. Journal of the observations. 


Band 

L. O. 

Date of 

Bandpass 

Flux 

Sky frequency 

On source 

Number of 

rms 


[GHz] 

observations 

calibrator 

calibrator 

coverage [GHz] 

time [min] 

antennas 

[mJy] 

B3 

92 

25-Apr-2012 

3C273 

Titan 

84-87.5, 96.2-99.7 

27 

19 

2 

B3 

95 

21-Nov-2012 

3C273 

3C273 

87.2-90.7,99.5-103 

17 

25 

2 

B3 

98 

21-Nov-2012 

3C273 

3C273 

90.5-94, 102.7-106.2 

10 

25 

2 

B3 

102 

25-Apr-2012 

3C273 

Titan 

93.7-97.2,106-109.5 

27 

19 

2 

B3 

105 

21-Nov-2012 

3C273 

3C273 

97-100.5,109.2-112.7 

18 

25 

2 

B6 

224 

31-Aug-2012 

3C279 

Titan 

214.3-217.8,229.3-232.8 

7 

23 

6 

B6 

253 

01-Jan-2013 

3C279 

Titan 

253.2-256.7,268.2-271.7 

7 

24 

6 

B6 

256 

01-Jan-2013 

3C279 

Titan 

246.7-250.2,261.7-265.2 

7 

24 

6 

B6 

259 

01-Jan-2013 

3C279 

Titan 

250-253.5,265-268.5 

7 

24 

6 

B7 

283 

2O-N0V-2OI2 

3C273 

Titan 

275-278.5, 287-290.5 

5 

25 

4 

B7 

286 

2O-N0V-2OI2 

3C273 

Titan 

278.2-281.7, 290.2-293.7 

5 

25 

4 
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Fig. 1. Total observed spectrum, extracted by fitting the visibilities with a point source. Each scan is plotted in a different color. The dashed line 
represents a 3 rd order fit to the line-free channels, and is used as continuum level in the analysis. 


it to be used for continuum removal in our fitting procedure. We 
choose instead to fit the line-free channels in the three bands with 
a polynomial. We obtain the best fit with a 3'^‘^ order polynomial, 
which is shown in Figures [T] and We will use this estimate in 
the following analysis of the molecular emission. 


By removing the estimated continuum flux from the spec¬ 
trum we can derive the contribution of molecular emission to 
the total observed flux in each band. We find that the molecular 
emission contributes for 15% of the total flux in band 3, and for 
27% in b ands 6 and 7. These values are very similar to the 28% 
found by iMartm et al.l (1201 ih for observations at 1 mm in the 
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Continuum at mm wavelengths 



Frequency [GHz] 


Fig. 2. Continuum flux density of NGC 4418 at mm wavelengths. 
The red diamonds represent the average value of line-free channels 
in the ALMA bands 3, 6, and 7. Sub-mm con tinuum levels from 
I Sakamoto et al.l (1201 Oli and I Sakamoto et al.l ll2013h are shown as blue 
triangles. Error bars were calculated as the quadrature sum of the spec¬ 
trum rms and a calibration uncertainty of 10%. The dot-dashed, and 
dashed lines represent the power-law fit of data points above 200 GFlz 
(a = 3 ± 0.6) and below 270 GFlz {a = 2 ± 0.4), respectively. The 
red solid line shows the best fit to the ALMA baseline with a 3^'*- 
order polynomial. The best fit parameters are A=2.7 x 10“’’ Jy/GFlz^, 
B=6.2 X 10 ** Jy/GHz^, C=3.2 x lO^^ Jy/GHz, D=2.3 x 10 ^ Jy. 


ULIRG Arp 220. These results clearly show that contamination 
by molecular emission, if not properly removed, could be a seri¬ 
ous issue for studies of continuum emission in the compact cores 
of (U)LIRGs. 

3.2. Line identification and fit 

The heavy blending of emission features makes the identification 
of individual molecular species a challenging task, especially in 
bands 6 and 7. In order to identify the emitting molecules we use 
a mix of local-thermal-equilibrium (LTE) and non-LTE (NLTE) 
fitting methods to fit the whole spectrum in the three ALMA 
bands. In the following paragraphs we describe the fitting routine 
and its results. 


3.2.1. Population (diagrams 

As a first step to the fit of the whole spectrum, we perform a pop¬ 
ulation diagram analysis for the molecules which have at least 
two non-blended lines in different ALMA bands. The integrated 
flux of each line was obtained by fitting a Gaussian profile to 
the spectrum. In the LTE approximation we assume the lines to 
be optically thin and the relative intensities to be described by a 
single ex citation temperature T r x- Lor each rotational transition, 
following iGoldsmith & Langeii(ll999l) . we can write the column 
density of the upper energy state as 


_ ^nkv^W 
hc^Aut 


(1) 


where v is the transition’s frequency, Au^ is the Einstein coeffi¬ 
cient, W is the integrated brightness temperature, h is Planck’s 
constant, k is Boltzmann’s constant, and c is the speed of light in 
vacuum. The integrated brightness temperature W is related to 
the observed integrated flux density S hne by 


lT = 5li„eX 


2kLl,’ 


( 2 ) 


with A the wavelength of the emission, and fis the solid angle 
subtended by the source. At LTE, we have 


N 

Nu = —gue 


-EJT„, 


( 3 ) 


where N is the total molecule’s column density, Z is the partition 
function, and gu and are respectively the degeneracy and the 
energy of the upper level. By fitting Eq.[3]to the derived Nu we 
can derive Tex and N for each molecule. Once we obtained Nu in 
the optically thin appr oximation, a correction for op acity effects 
was applied, following iGoldsmith & Langeil(ll999l) : 




3h Av Z 



^thicls _ ^thin 

ly u 4 V y /N 


1 


( 4 ) 


where ^ is the dipole moment, Av the full width at half maximum 
(LWHM) of the emission line, and Bq the rotational constant of 
the molecule. 

The transition parameters (v, Eu, Aut, gu) were taken from 
the SplataloguJ^ database, and were mainly deri ved from the 
Colog ne Database for Molecular Spe ctroscop'^ ([M uller et ^ 
I 2 OO 5 I) and Jet Propulsion Laborator^ (iPickett et alJI 19981) cata- 
logs. The partition function was calculated by summation on all 
the available transitions as Z = l,ugue^^'‘^^^°'- 

We take as source size of the emission the upper limit 
of 0".4 to the HCN and HNC J=3-2 emission observed by 
I Sakamoto et al.l (1201 3l) and ICostagliola et^ (l2013l) . Notice that 
this value is compatible with the upper limit to the source size of 
HNC 1-0 in our data. 

We fit transitions from twelve different species; HCN, HNC, 
HCO+, CS, C^^s, i3cs, SiO, HC 3 N, vibrationally excited HC 3 N 
in the v7=l, v 6 =l, and \l-2 states, and HC 5 N. The results are 
summarized in Pig. [3] and Table |2| We find excitation tempera¬ 
tures ranging from 7 to 140 K, and column densities from 10'^ to 
10 '^ cm“^. The scatter in excitation temperatures is larger than 
the uncertainties and may be due to a steep temperature gradient 
or NLTE excitation effects. 

Standard high-density tracers as CS, HCN, HNC, and HCO^ 
have low excitation temperatures (<20 K) and high opacities 
(>1), while HC 3 N transitions have moderate to low opacities and 
higher excitation temperatures, of the order of 100 K. Given the 
uncertainty on the source size and beam filling factor, the de¬ 
rived column d ensities have to be conside red as lower limits to 
the real value (iGoldsmith & LangeJll999l) . This is particularly 
true for molecules with moderate to high opacities. 

3 www.spiataiogue.net 
www.astro.uni-koeln.de/cdms 
^ spec.jpl.nasa.gov 
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Table 2. Gaussian fitting of blend-free lines and results of the population diagram analysis. 


Molecule 

Transition 

Erequency 

[GHz] 

0(1) 

^line 

[Jy kms^'] 

[kms-'l 

’J-'O) 

-'b 

[K] 

^ ex 

[K] 

[cm“^] 

,^(6) 

HCN 

J=l-0 

88.63 

9+1 

140+10 

65+7 

8+2 

2.2+0.5x10*’ 

>10 


1=3-2 

265.89 

50+5 

150+10 

35+4 



>10 

HNC 

1 = 1-0 

90.66 

5.2+0.7 

110+15 

45+6 

10+3 

1.7+0.5x10**’ 

5.9 


1=3-2 

271.98 

38.5+0.3 

130+10 

30+2 



6.8 

HCO+ 

1 = 1-0 

89.19 

6.2+0.4 

110+7 

55+3 

6+1 

1 . 1 + 0 . 2 x 10 *® 

>10 


1=3-2 

267.56 

23.6+2 

150+12 

17+2 



6.6 

CS 

1 = 2-1 

97.98 

5+0.5 

116+12 

35+4 

15+3 

3.6+0.8x10*® 

2.3 


1=6-5 

293.91 

32+4 

150+16 

19+2 



1.7 


1 = 2-1 

96.41 

1+0.3 

100+30 

8+3 

16+10 

5+3x10*® 

0.6 


1=6-5 

289.21 

10+1 

115+15 

8+1 



0.8 

'^cs 

1 = 2-1 

92.49 

0.3+0.2 

90+10 

3+2 

22+7 

3+1x10*® 

0.1 


1=5-4 

231.22 

L6+0.7 

90+40 

3+1 



0.1 


1=6-5 

277.45 

3.3+0.8 

150+40 

2+0.5 



0.1 

SiO 

1 = 2-1 

86.85 

L2+0.3 

90+20 

13+3 

9+5 

3+1x10*® 

1.5 


1=5-4 

217.10 

4+1 

130+30 

5+1 



0.8 

HC 3 N 

1=10-9 

90.98 

4.0+0.8 

110+20 

34+7 

65+4 

1 + 0 . 1 x 10 *’ 

0.5 


1 = 11-10 

100.08 

4.6 +0.6 

110+15 

31+4 



0.4 


1 = 12-11 

109.17 

5.1 +0.8 

120+20 

27+4 



0.4 


1=24-23 

218.32 

15.2 +2.5 

150+25 

16+3 



0.2 


1=28-27 

254.70 

23.8 +2.2 

150+15 

18+2 



0.3 


1=29-28 

263.79 

17.4 +2.2 

135+15 

14+2 



0.2 


1=30-29 

272.88 

16.2+1.9 

135+15 

12+1 



0.2 


1=31-30 

281.98 

19.4 +2.0 

150+15 

12+1 



0.2 


1=32-31 

291.07 

19.4+1.7 

150+15 

12+1 



0.2 

HC3N,v7=1 

J=10-9,l=lf 

91.33 

0 . 8 + 0.2 

90+25 

8+2 

98+6 

1.4+0.1x10*® 

0.1 


J=ll-10,l=lf 

100.47 

1.3+0.3 

105+30 

9+2 



0.1 


l= 12 -ll,l=lf 

109.60 

L2+0.3 

90+30 

8+2 



0.1 


1=24-23,1= If 

219.17 

7+ 1.5 

130+30 

9+2 



0.1 


1=28-27,1= If 

255.70 

10+2 

150+25 

8+1 



0.1 


1=29-28,1= If 

264.82 

9+2 

130+30 

7+2 



0.1 


J=31-30,l=lf 

283.07 

11 + 1 

145+10 

7+1 



0.1 


l=32-31,l=lf 

292.20 

10+2 

120+20 

7+1 



0.1 

HC3N,v6=1 

J=10-9,l=le 

91.13 

0.50+0.4 

140+100 

3+2 

99+21 

8 + 2 x 10 *’ 

«1 


J=ll-10,l=le 

100.24 

0.4+0.3 

140 + 100 

2+1 



«1 


l= 12 -ll,l=le 

109.35 

0 . 6 + 0.2 

120 + 40 

3+1 



«1 


1=29-28,1= le 

264.22 

5+1 

140 + 30 

4+1 



«1 

HC3N,v7=2 

J=10-9,l=2e 

91.56 

0.4+0.2 

70+40 

5+3 

139+27 

1.9+0.4x10*’ 

«1 


J=ll-10,l=2e 

100.71 

0.9+ 0.3 

105+30 

7+2 



«1 


J=12-ll,l=2e 

109.87 

1.1+0.3 

120+35 

6+2 



«1 


l=28-27,l=2e 

256.31 

13.1+2.5 

170+30 

8+2 



0.1 


l=32-31,l=2e 

292.91 

11 . 2 + 2.6 

170+40 

6+1 



«1 

HC 5 N 

1=32-31 

85.20 

0.6+0.5 

150+100 

5+4 

36+22 

1.2+0.7x10*® 

0.1 


1=33-32 

87.86 

0 . 6 + 0.2 

130 +45 

4+1 



0.1 


1=34-33 

90.52 

0 . 2 + 0.1 

100 +50 

2+1 



«1 


1=37-36 

98.51 

0.6+0.3 

130+65 

4+2 



0.1 


1=38-37 

101.17 

0.4+0.1 

90+30 

4+1 



0.1 


1=39-38 

103.83 

0 . 2 + 0.1 

90+50 

2+1 



«1 


1=42-41 

111.82 

0.4+0.1 

90+35 

3+1 



0.1 


(1) Flux density integrated over the fitted Gaussian profile; (2) Gaussian line width; (3) Brightness temperature assuming a source size of 0" .4. (4) Excitation 
temperature from population diagram fit; (5) Column density from the population diagram fit; (6) Opacity of the transition. 


3.2.2. LTE/NLTE fit of the whole spectrum 

In order to identify the fainter, heavy blended emission, we fit the 
whole spectrum with a combined LTE/NLTE mode l, in a fashion 
similar to what was done for the ULIRG Arp 220 bv lMartln et al.l 


(1201 ll) . We include in the fit all the molecules ident ified so far in 
theextragalactic interstellar medium, as reported in iMartm et al.l 
(1201 ll Table 5), for a total of 46 species and 23 isotopologues. 
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Fig. 3. Population diagrams for molecules with non-blended emission lines. The diamonds represent the values of Nu/gu for each transition, 
derived applying the opacity correction. The solid line shows the fit to the opacity-corrected data, while the dashed line shows the fit assuming 
optically thin emission. The fitted values for column density (N) and excitation temperature (Tex) are reported, together with the maximum opacity 
of the observed transitions. The error bars show 3-cr uncertainties. 


The fitting procedure is composed of the following steps: 7j 
a synthetic LTE spectrum is produced for each molecule, assum¬ 
ing a fiducial column density of 10^^ cm“^, and an excitation 
tempe rature of 80 K (based on observations bv ICostagliola et all 
1201 3h : 2) the molecules are ordered from the brightest to the 
faintest based on the flux density from the synthetic spectrum in 
bands 6 and 7; a synthetic spectrum for the first molecule is 
fitted to the continuum-removed observed spectrum by mini¬ 
mization; 4) The best fit spectrum is removed from the data and 
the fit parameters stored. Steps 3 and 4 are repeated iteratively 
for all the molecules in the list. 


For the molecules with collisional coefficients available in 
the LAMBDA database (CS, S O, SiO, HCN HNC , HCO+, 
N2H+, H2CO, CH3CN, CH3OH, ISchoier et al.ll200-5h . both an 
LTE and an NLTE synthetic spectrum was fit to the data, while 
for the other molecules only an LTE fit was possible. Eor the 
NLTE fit, a grid of 100,000 synthetic s pectra was produced wit h 
the RADEX radiative transfer code (Ivan der Tak et alJ [2007h . 
with molecular hydrogen densities ranging from 10^ to lO^cm^^, 
kinetic temperatures between 5 and 300 K, and molecular col¬ 
umn densities between 10*^ and lO'^ cm“^. The range in ki¬ 
netic temperatures was limited by the available collisional co¬ 
efficients, which in most cases have been computed up to 300 K. 
The LTE synthet ic spectra were prod uced following the proce¬ 
dure described in iMartm et al.l (1201 ll) . 


We considered a fixed line width of 120 km s', which is the 
average of the values found by Gaussian fitting the blend-free 
lines (see Section [3.2.1l and Table |2]). Before calculating the^^, 
the synthetic spectra were red-shifted assuming a single LSR ve¬ 


locity for NGC 4418 of 2100 kms“*. We assumed all the emis¬ 
sion to be emerging from a region of 0".4 in diameter, which is 
the upper limit to the size of mole cular core of NGC 4418, as ob - 
served by ISakamoto et al.l (1201 3h and ICostagliola et alJ (l2013l) . 
In both the LTE and NLTE fit, we use the 3’^“' order fit of the 
line-free channels as continuum level, as explained in Section 

ITT] 


3.3. Results from the line identification and fit 

A list of the identified molecular species is reported in Table [3] 
We define as detected those molecules which have more than 
one transition included in our spectral scan and are well fitted by 
our LTE or NLTE models. An exception to this rule are bright 
lines from well-known ISM species or their isotopic variants, 
such as *^CO, H'^CN, HN'^C, and H'^CO^, for which one line 
is enough for a successful identification. We define as tentative 
detections those molecules with only one detected transition in 
the scan (e.g., HCN,v2=l), or those for which no satisfactory fit 
could be obtained (e.g., CH3OH). Tentative detections are shown 
in italics in Tabled 

We identify 317 emission lines above the three sigma level 
from a total of 45 molecular species, including 15 isotopic sub¬ 
stitutions and six vibrationally excited variants. The emission 
lines above three sigma which were fitted by our model are listed 
in Tables lATI to lA.hl The flux recovered by the fit is 92 % of the 
total molecular emission, resulting in a normalized;\^^ of 5 for the 
whole spectrum. This somewhat high value of the is mainly 
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Table 3. Summary of species detected in NGC 4418. Tentative detections are shown in italics. 


2 atoms 3 atoms 

4 atoms 

5 atoms 

6 atoms 

7 atoms 

CS 

HCN \ 

P-H 2 CO 

HC 3 N ) 

CH 3 CN 

CHiCCH 

'^cs 

H'^CN / 

0 -H 2 CO 

HCC'^CN i 

CH 2 OH 

HC 5 N 


HCN,v2=l 

c-HCCCH 

H'^CCCNj 




HNC 1 

H 2 CS 1 

HC3N,v6=1 



'^co 

1 HN'^C j 

//>^CS J 

HC3N,v7=1 




1 HNC,v2=l 


HC3N,v6=1,v7 

= 1 


CN 

HCO+ ) 


HC3N,v7=2 



NS 

Hi^CO^ i 


CH 2 NH 



SO 

j 


NH 2 CN 



3450 








CCH 

HCS+ 

CCS 

N2H+ 


Vibrationally excited HC^N 


15r 


14 


v 6 =l 


13 


v =0 


-s? n 

Z 

^ 11 
10 - 
9 


v7=l 


v7=2 


^_0 J=10-9, T^^=330 (80) K 
o o o J=11-10, T^^=340 (80) K 
A J=12-11,T =400(100)K 


200 


400 

E.. [K] 


600 


800 


Fig. 4. Population diagram for the vibrational excitation of HC 3 N. For 
the three rotational levels J=10 (diamonds), J=ll (circles), J=12 (tri¬ 
angles) we compare the relative populations of the v=0, v7=l, v7=2, 
and v 6 =l vibrationally excited states. The vibrational temperature fit is 
shown as a dashed, dotted and dash-dotted line, for the J=10-9, 11- 
10 , and 12-11 transition respectively. The three fits result in similar 
vibrational temperatures in the range 330-400 K. 


due to unidentified lines and systematic uncertainties, which are 
discussed in Section[33] 

The results of the fit are summarized in Table 01 The fit de¬ 
tails for individual molecules are reported in SectionlAl 


3.4. Limitations of the fit and main uncertainties 

The confidence intervals reported in Tables |2] and 0] represent 
one-sigma statistical uncertainties derived from thex^ minimiza¬ 


tion. Other systematic uncertainties have to be considered for a 
correct interpretation of the fit results. 

Source size and beam filling factor: Because the emission 
is unresolved in the ALMA beam, the main uncertainty in our 
LTE/NLTE fit is the source size. In our conversion from flux 
density to brightness temperature we assume a source size of 
0."4, which is the upper limit to the diameter of the region 
containing the bulk of the molecular emission in NGC 4418 
(ISakamoto et all l2()T^ ICostagliola et d1 l2013h . If we assume 
that all the molecular emission is coming from this compact re¬ 
gion, the derived brightness temperature is a lower limit to the 
true value. At LTE this translates into underestimated column 
densities or excitation temperatures, respectively for optically 
thin or optically thick transitions. 

Our fit also assumes that all the molecular species have 
the same beam filling factor, which is unlikely, given the steep 
temperature and density gradi ents expected in this galaxy (e.g., 
iGonzalez-Alfonso et alJl2()T^ . The chemical properties derived 
for the ISM of NGC 4418 have thus to be intended as an average 
of the different environments contained in the ALMA beam. The 
same limitation is true for most of the extragalactic spectral scans 
observed so far, and will be possible to overcome only with high- 
resolution observations. Given its extreme compactness, high an¬ 
gular resolution observations (e.g. with the full ALMA array in 
an extended configuration) will be required to further resolve the 
compact mm-line emission region. 

Line blending: Severe line confusion in bands 6 and 7 makes 
the fit non-univocal for the species with no strong blend-free 
lines detected. This problem is somewhat mitigated by the large 
frequency coverage of our scan, but its effects are evident in 
some heavily blended regions as a poor fit to the data. Based on 
the relative intensity of the blended lines, we estimate the con¬ 
tribution to the uncertainty on the fitted parameters to be of the 
order of 20 %. 

Opacity: Eor optically thick transitions the intensity of the 
lines depends only marginally on the column density. Eor the 
most abundant molecules, such as CS, HCN, HNC, and HCO^, 
which have opacities greater than one in most of the observed 
transitions, the column density is thus poorly determined. Eor 
this reason we choose not to derive isotopic ratios for optically 
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Table 4. Results from the LTE and NLTE fit of the whole spectrum. 



LTE 

NLTE 

Molecule 

Tex [K]"> 

N [cm-2]m 

n(H 2 ) [cm-®]® 

Tkin [K]® 

Tex®-'”” [K]<=> 

N [cm-®]® 

CS 

20+1 

6 + 0.6 xl 0 ‘® 

10 ® 

55+2 

47-30 

1 + 0 . 2 x 10 '® 


47 

L5+0.5 xlO'^ 

- 

- 

- 

- 


47 

3+1 xlO's 

- 

- 

- 

- 


47 

8 + 2 x 10 '^ 

- 

- 

- 

- 


70 

5+2x10'® 

- 

- 

- 

- 


70 

6+1 xlO'^ 

- 

- 

- 

- 

CN 

70 

5+0.5 xlO'^ 

- 

- 

- 

- 

NS 

350+150 

8+4x10"’ 

- 

- 

- 

- 

SO 

>20 

lO'^-lO'® 

10 ®- 10 ® 

50-200 

40-100 

2 - 20 x 10 '® 


>20 

10'5-ioi6 

- 

- 

- 

- 

SiO 

10+5 

3+1 xlO'^ 

43 

0 

V 

>20 

22-10 

V 

0 

2 ‘’SiO 

22 

2+1 xlO's 

- 

- 

- 

- 

30SiO 

22 

2+1 xlO'^ 

- 

- 

- 

- 

HCN 

7+1 

1.5+0.1 xlO'5 

< 10 ® 

>100 

23-15 

> 10 '® 

h^^cn 

65 

1+0.5 xlO'® 

- 

- 

- 

- 

HCN,v2=l 

350 

3+1 xlO'^ 

- 

- 

- 

- 

HNC 

8+1 

8+0.5 xlO'5 

43 

0 

V 

>100 

45-40 

V 

0 

HN'^C 

45 

1+0.5 xlO'5 

- 

- 

- 

- 

HNC,v2=l 

450 

1.5+0.5 xlO'® 

- 

- 

- 

- 

HCO+ 

6+1 

6+0.5 xlO'5 

< 6 x 10 ® 

> 100 

80-50 

> 10 '® 

h'3co+ 

80 

2.5+1 xlO'^ 

- 

- 

- 

- 


<20 

>10'4 

- 

- 

- 

- 

H 2 S 

70 

1+0.5 xlO'^ 

- 

- 

- 

- 

CCH 

70 

2+1 xlO'^ 

- 

- 

- 

- 

HCS+ 

20+10 

8+4x10'^ 

- 

- 

- 

- 

CCS 

20+10 

L5+0.5 xlO'® 

- 

- 

- 

- 

N 2 H+ 

30+5 

5+1 xlO's 

10 ® 

23+2 

23-23 

5+2 XlO'® 

P-H 2 CO 

350+100 

5+2x10'® 

10 ® 

>300 

20-100 

7+2x10'® 

0 -H 2 CO 

350+100 

8+4x10'® 

10 ® 

>300 

100-500 

1 + 0 . 2 x 10 '® 

c-HCCCH 

180+100 

1+0.5 xlO'^ 

- 

- 

- 

- 

H 2 CS 

35+10 

1+0.5 xlO'® 

- 

- 

- 

- 

H^^CS 

35 

5+2x10'® 

- 

- 

- 

- 

HC 3 N 

70+5 

1.8+2 xlO'® 

- 

- 

- 

- 

H^^CCCN 

70 

2+1 XlO'® 

- 

- 

- 

- 

HCC^^CN 

70 

1+0.5 XlO'® 

- 

- 

- 

- 

HC3 N,v6=1 

100+50 

4+3 XlO'® 

- 

- 

- 

- 

HC3N,v7=1 

100+25 

1+0.2 XlO'® 

- 

- 

- 

- 

HC3N,v6=1,v7=1 

70+30 

3+2x10'® 

- 

- 

- 

- 

HC3N,v7=2 

100+33 

8 + 2 x 10 '® 

- 

- 

- 

- 

CH 2 NH 

45+30 

L5+0.5 xlO'^ 

- 

- 

- 

- 

NH 2 CN 

30+20 

5+2 XlO'® 

- 

- 

- 

- 

CH 3 CN 

>1000 

8+3 xlO'® 

10 "' 

>300 

10-20 

1.5+0.3x10'® 

CHyOH 

60 

< 1 x 10 '® 

10 ® 

60+10 

20-40 

2 + 1 x 10 '® 

CHiCCH 

350+100 

1+0.5 xlO'^ 

- 

- 

- 

- 

HC 5 N 

65+7 

1+0.4x10'® 

- 

- 

- 

- 


(I) Excitation temperature from LTE fit; (2) Column density from LTE fit; (3) Molecular hydrogen density from NLTE fit; (4) Kinetic temperature from NLTE fit; (5) 
Maximum excitation temperature at 3 and 1 mm from NLTE fit; (6) Column density from NLTE fit. 


thick species, and we will not include these molecules in the 
following abundance analysis. 

NLTE effects: An NLTE analysis of the emission was only 
possible for the molecules included in the LAMBDA database, 
while for the others only an LTE ht was performed. The NLTE fit 


results in a wide range of molecular hydrogen densities, it is thus 
likely that NLTE effects may affect many of the other species. 
The LTE assumption implies densities greater than the critical 
density for the transitions, which in general increases with the ro¬ 
tational quantum number. At lower densities, the de-population 
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of the high-energy levels is interpreted by our LTE code as an 
underestimated kinetic temperature. This results evident when 
comparing the LTE and NLTE fits for HCN, HNC, and HCO+. A 
more precise analysis will be possible only when the collisional 
coefficients for all the relevant species will become available. 

Line width: In order to fit the heavily blended lines in bands 
6 and 7, we assume a single line width of 120 kms^*, which is 
an average value derived from the blend-free transitions of Table 
13 However, the measured line widths have a significant scatter, 
with a trend of an increasing EWHM with critical density (see 
Fig. IS. The measured line widths have a standard deviation of 
17 km s ', which translates into an error on the fitted integrated 
intensities of the order of 15%. 

Most of the systematic uncertainties in our fitting strategy re¬ 
sult in a potential underestimation of molecular column density 
or kinetic temperature. The effect of other uncertainties (such as 
those introduced by line blending or variations in line width) on 
the fit results does not exceed 20%, comparable to the nominal 
accuracy of the data calibration. 

4. Discussion 

4.1. A vibrationally excited spectrum 

The spectrum of NGC 4418 appears to be dominated by bright 
emission from vibrationally excited HC3N, which was detected 
in the v7=l, \l-2, v6=l, and (v7=l,v6=l) variants. The detec¬ 
tion of emission from the (v7—l,v6—l) state of HC^N is the first 
obtained in an extragalactic object. We also detect the J=3-2 
transitions of vibrationally excited v2=l HCN and HNC. Emis¬ 
sion from vibrationally excited HC3N, HCN, and HNC amounts 
to 19% of the total flux recovered by our fit. 

The excitation energies of the HC3N vibrational levels range 
from 300 to 900 K, corresponding to transition wavelengths in 
the far infrared (EIR), from 16 to 48 pm. The population analysis 
of vibrationally excited HC3N shown in Eig.|4]results in a vibra¬ 
tional temperature of 330-400 K. Even higher values of 350 and 
450 K are found respectively for HCN and HNC when compar¬ 
ing the intensities of the 1=3-2 transitions of the v=0 and v2=l 
vibrational states. 

The critical densities of the vibrational EIR transitions of 
HC3N, HCN, and HNC are high, greater than 10^ cm“^. We 
can thus assume that, for typical ISM densities (10^-10^ cm“^), 
radiative excitation dominates the population of the vibrational 
levels of HC3N. Under this assumption, the relative population 
of the vibrationally excited levels is a direct probe of the EIR ra- 
diation field in the core of NGC 4418 (e.g., ICostagliola & Aalt^ 
I2OIOI) . 

In the limit of perfect coupling between the EIR radiation 
field and the vibrational levels of the three molecules, the derived 
vibrational temperature is a good approximation to the tempera¬ 
ture of the infrared source in the core of NGC 4418. When allow¬ 
ing for a non-uniform coverage of the IR source, inefficient cou¬ 
pling, or contamination from collisional excitation, the derived 
value has to be interpreted as a lower limit to the temperature of 
the exciting IR radiation. 

The only extragalactic object where a similar vibrationally 
excite d spectrum has bee n detected is the ULIRG Arp 220, 
where iMartfn et alJ (1201 ih find remarkably similar vibrational 
temperatures (300^00 K) for HC3N. Both objects have strong 
IR fluxes, very compact molecular cores (<20 pc), and high 
molecular columns (N(H2)>10^"^ cm“^); a vibrationally excited 
spectrum may thus be the signature of compact, deeply buried IR 
sources in extragalactic objects. 


4.2. Temperature and density structure 

Our fit of the ALMA spectral scan finds kinetic temperatures 
ranging between 20 and 350 K, and molecular hydrogen densi¬ 
ties between 10^ and 10^ cm“^. 

The vibrational excitation of HC3N, HCN, and HNC, re¬ 
quires the presence of a deeply buried, optically thick source, of 
brightness temperature higher than 350 K. If we assume all the 
IR flux of NGC 4418 to be coming from such a source, this must 
be extremely compact, with a size smaller than 5 pc. A com¬ 
pact IR sourc e of similar temperature and si ze is also required 
by models by iGonzalez- Alfonso et al.l (1201 2h in order to fit the 
mid-IR continuum. 

The single sphere geometry is a simplifying assumption and 
our data do not distinguish a single 5 pc sphere from other ge¬ 
ometries such as a disk or a group of smaller patches distributed 
on a larger scale, as long as they have the same total area. We 
note that a radio source of ~5 pc in diameter has a ctually been 
detected by EVN observations (IVarenius et al.ll2014l) at less than 
10 milli arcseconds fr om the peak of the 860 pm continuum 
(ISakamoto et al.ll2013h . This source may be responsible for the 
compact IR and vibrationally excited emission, however this sce¬ 
nario needs spectroscopic observations at higher spatial resolu¬ 
tion in order to be confirmed. 

Very high rotational temperatures, exceeding 300 K, are 
found for molecules such as NS, H2CO, CH3CN, and CH3CCH. 
The observed transitions have critical densities lower than 
10® cm“^ and we can assume the excitation to be dominated 
by collisions under normal dense-ISM conditions. This is the 
first time such high rotational temperatures have been found in 
NGC 4418, mainly thanks to the wide frequency coverage of our 
scan. These measurements reveal a warm, >300 K gas compo¬ 
nent, which may be associated with the optically thick IR source 
producing the vibrationally excited emission. 

The ISM in the core of NCG 4418 appears to be a multi¬ 
phase environment, where we can identify three main temper¬ 
ature components; one cold component at 20-50 K, traced by 
CS, SiO, HCS+, CCS, N2H+, H2CS, CH2NH, and NH2CN; 
one warm component at 60-100 K, traced by CO, CN, HCN, 
HNC, HCO+, H2S, CCH, HC3N, CH3OH, and HC5N; and one 
hot component at >300 K, traced by the vibrationally excited 
HC3N, HNC, and HCN, and by collisionally excited NS, H2CO, 
CH3CN, and CH3CCH. 

Our analysis is consi stent with the resu lts _ by 

IGonzalez- Alfonso et al.l (1201 2l) and ICostagliola et al.l (1201 3ll . 
who modeled Herschel/PACS and SMA observations with a 
layered density and temperature structure. 

Additional evidence of a radial density and temperature gra¬ 
dient comes from the analysis of the velocity dispersion of the 
blend-free lines of Table |3 Eigure |5] shows that the line width 
increases with critical density and upper state energy of the tran¬ 
sitions. If we assume the velocity dispersion to be radially in¬ 
creasing towards the center of the galaxy, as it would be the 
case for a Keplerian rotation, we can interpret these trends as 
a radial gradient of density and temperature. This result is in 
agreement with our previous S MA observations (ISakamoto et al.l 
l2013t l^stagliola et al.ll2013h . which found a similar trend for 
the molecular emission at 1 mm and 850 pm. These studies also 
revealed that emission from transitions with high critical densi¬ 
ties tend to be more centrally concentrated, supporting the exis¬ 
tence of a radial density gradient. 

Most of the transitions considered in Eig. |5] however have 
opacities greater than one and we cannot exclude the observed 
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Fig. 5. Full width at half maximum versus critical density and upper state energy. The line width was derived by a Gaussian fit of the blend-free 
lines (see Table |2](. The critical density (nc) was calculated from the Einstein’s and collisional coefficients in the LAMBDA database. The upper 
state energies (E^,) of the transitions were taken from Splatalogue. 


line broadening to be caused by radiative trapping. To confirm 
this scenario, we will need spatially-resolved observations. 

4.3. Extragalactic HC 5 N 

Emission from HC 5 N has been observe d towards several Galac¬ 
tic hot cores (e.g., iGreen et^ l2014l) . but until recently the 
molecule had not been detected outside the Galaxy. Aladro et 
al. {in prep.) report the detection of eight HC 5 N transitions in 
the starburst galaxy NGC 253, to our knowledge NGC 4418 is 
the only other extragalactic object where the molecule has been 
detected. 

From the eleven detected transitions we derive an excita¬ 
tion temperature of 60-70 K, and a column density of 0.5- 
2x10'® cm“^. Assuming a column density of H? towards th e nu¬ 
cleus of NGC 4418 of 10^4 cm^^ dCostagliola et akllMl . this 
corresponds to an abundance of 0.5-2x10^*. This value is about 
one order of magnitude hig her than what fo u nd in N GC 253. 

Chemical models by IChanman et al.l (l2009ll show that 
cyanopolyynes can reach high abundances in hot cores, thanks 
to gas-phase neutra-neutral reactions. The main formation path 
of HC 5 N for gas temperatures of the order of 100 K is 

C 2 H + C 2 H 2 H 2 CCCC + H 

CN + H 2 CCCC HCsN + H, (5) 

which for molecular hydrogen densities close to 10 ® cm ® can 
lead to detectable column densities (N>10*"' cm“^). 

The abundance of HC 5 N is directly linked to acetylene, 
which is highly abundant in the earlier stages of the core evo¬ 
lution when it is evaporated from grain mantles. The enhance¬ 
ment of HC 5 N in the gas phase is extremely short lived, it abun¬ 
dance increasing and decre asing over a pe r iod of the order of 
10 ^- 10 ® yr. For this reason. IChapman et al.l (l2009t) suggest that 
the molecule could be used as a chemical clock in hot-cores. 
The fact that HC 5 N has been detected so far only in two exter¬ 
nal galaxies may be due to its short lifetime in the gas phase. 
The detection of HC 5 N in NGC 4418 may thus be coming from 
an extremely young, d ust- embedded starbur s t. Obs ervations by 
I Sakamoto et al.l (l2013h and ICostagliola et^ (|20 1 3h indeed sug¬ 
gest that if a starburst is powering the central 100 pc of the 


galaxy, this should be younger than 10 Myr. The fact that the 
molecule has been detected in NGC 253 but not in M 82 supports 
the interpretation that its abundance may be a good indicator of 
the starburst evolution. 


4.4. Standard line ratios: HCN, HNC, and HCO^ 


Line ratios of the dense gas tracers HCN, HNC, and HCO'^ 
have often been used in the literature to investigate the nature 
of the dominant energy source in LlRGs, either a compact star- 
burst or an AGN (e.g., Imanishiet^ 20071 : lAalto et ^12007 at 
iKrips et ^l2008HCostagliola et al.ll20lT ). 

HCO^/HCN ratios: The HCO^/HCN J=l-0 line ratio has 
been proposed as a diagnostic tool to d istinguish between 
AGN- and starburst-powered galaxies (e.g., iKohno et al.ll2001t 
]2007h . Interferometric observations by 
reveal that Seyfert galaxies have lower 


Imanishi et al. 


Kohno et al.l ( 200111 


200 


HCO^/HCN ratios compared with starbursts. Seyfert galaxies 
with starburst-like HCO^/HCN ra tios are generally int erpreted 
as mixed A GN-starburst object s (llmanishi et al.ll2007h. S ingle 
dish surveys iKrins et ^ (l2008h : ICostagliola et al.l (1201 3h find 
that starbursts and LIRGs have, on average, higher HCO^/HCN 
ratios than AGNs, which in this scenario would suggest a dom¬ 
inant UV-dominated (PDR) chemistry. The theoretical explana¬ 
tion of the observed low HCC )^/HCN in AGNs i s how ever still 
being debated. Early work bv iLenp & Dalgarn^ (Il996h suggest 
that the observed trend would be due to an enhancement of HCN 
abundance in the X-ray-dominated region surrounding an AGN. 
However, more recent chemical models find that the HCO'^/HCN 
abundance ratio is a strong functio n of density and ionizing radi¬ 
ation, (e.g., iMeiierink et al.l2d0^ . and that a low line ratio could 
be either caused by low-density(<10® cm XDRs or dense 


(> 10 ® 


PDRs. Also, given the different critical densities 


of the J=l-0 transition of HCN and HCO^, localized variations 
in the ratios could be explained by variations in density (e.g., 
iMeier & Turneijl2012h . 

In our spectral scan of NGC 4418 the HCO^/HCN line ratio 
varies between 0.8 to 0.5 for the J=l-0 and J=3-2 transitions re¬ 
spectively. The abundance ratio can be estimated by H'®CN and 
H'®CO+ J=l-0 transitions. Assuming optically thin emission. 
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and an excitation temperature of 7 and 6 K, we find LTE column 


densities of 2+1x10*'’cm-^ and7±4xl0*^cm-^ forH'^CNand 
respectively. If we assume an extrag alactic ra¬ 

tio between 40 and 100 (IHenkel et al.ir2014l) . these correspond 
to column densities of (0.8-2)xl0*^ cm“^ for HCN and (2.8- 
7)xl0'*’ cm-2 for HCO+. The derived HCOVHCN abundance 
ratio in NGC 4418 is thus 0.35+0.2. 

The molecular gas density cannot be higher than 
10^ cm^^ because at higher densities the HCO*^ 3-2 tran¬ 
sition would be thermalized, while we observe a J=3-2/l-0 
brightness temperature ratio of 0.3. For such densities, the 
observed HCO*^/HC N line and abundance ratios are consistent 
with PDR models bv iMeiierink et al.l (l2007h . while XDR models 
would result in much higher ratios (between 3 and 5). 

However, high-resol u tion A LMA observations of NGC 1068 
by iGarcfa-Burillo et all ( |2014 hnd a HCO+(4-3)/HCN(4-3) 
brightness temperature ratio of ~0.4 in the region within 140 pc 
from the AGN, and a ratio of ~0.8 at the AGN position. These 
values are very close to the ones observed in NGC 4418 for 
the J=l-0 and J=3-2 transitions and provide a strong observa¬ 
tional evidence for the low HCO*^/HCN ratio to be associated 
with AGN activity. The higher ratio measured directly towards 
the AGN in NGC 1068 suggests that the relative intensity for 
the HCO*^ and HCN emission may be regulated by mechanical 
heating or shock chemistry from the AGN feedback rather than 
by X-ray dominated processes. 

HNC/HCN ratios: In the Galaxy, the HNC/HCN ob served 
line ratio ranges from 1/100 in hot cor es dSc hilke et ah 1992h 
to values as high as 4 in dark clouds (iHirota et al. 1998).The 


abundance of th e HNC molecule decreases with increasing 
gas temperature. iHirota et'aP (Il998h suggest that this may be 
due to the temperature dependence of neutral reactions, which, 
for temperatures exceeding 24 K, selectively destroy HNC in 
favour of HCN. Bright HNC e mission is commo n ly observed 
in ex t ragalactic objects (e.g .. lAa lto et al J 120021: IWang et akl 
l2004t iMeier & Turned [2005t iPerez-Beaunuits et al. 2007). In 
particular, over-luminous HNC J-3-2 is f ound in LIRGs (e.g., 
lAalto et ^l2007bl : ICostagliola et al.ir201 Ih . where gas tempera¬ 
tures, derived by the IR dust continuum and mid-IR molecular 
absorption, are usually higher than 50 K and can reach values as 
high as a few 100 K. Here, ion-molecule chemistry in PDRs may 
be responsible for the observed ratios. iMeiierink et al.l (l2007l) 
hnd that the HNC/HCN J=l-0 line ratio is enhanced in PDRs 
and can reach a maximum value of one for H 2 column densities 
exceeding 10 ^^ cm“^. 

In NGC 4418 we hnd HNC/HCN line ratios of 0.7 and 
unity for the J=l-0 and J=3-2 transitions, respectively. These 


values are in agreement with a pi 
tic survey by ICostagliola et alJ (y 


irevio us single-dish extragalac- 
201 If) , which hnds and average 
HNC/HCN J=l-0 of 0.5, with a maximum value of one, for a 
sample of 23 galaxies, including LIRGs. The observed ratios 
can be easily explained by P DR chemistry (IMeiierink & SpaaTii 
l2005UMeiierink et alJl2007l) . 

Although the line ratios of HCN, HNC and HCO*^ ob¬ 
served in NGC 4418 are compatible with PDR models, we 
advise extreme caution when interpreting these results. Ex¬ 
tragalactic surveys have shown that these line ratios show 
only little variations, of the order of a few, across different 
galaxy types and t h eir interpretation is still controversial (e.g. . 


Kohno et^ 20011: ^Ito et’aD I2007R IMeiierink et alJ 2007 : 
CostaglioHey^ joiil). Recent stud ies (e.g.. IMeiierink et alJ 


201 It iKazandiian et al.l 1201 2l 120151) show that including me¬ 


chanical heating (e.g., from star-formation or AGN feedback) 
and cosmic rays (e.g., from supernovae) in a PDR network can 


change the line ratios of HCN, HNC and HCO*^ by more than 
one order of magnitude. The observed values may be thus due to 
a combination of factors which are very difficult to disentangle 
in spatially unresolved observations. 

4.5. Chemical composition of NGC 4418 

In Fig | 6 ] we compare the molecular abundances found for 
NGC 4418 with those observed in other Galactic and extra¬ 
galactic objects. We choose to include in the analysis only 
molecules with low to moderate opacities, and for which reli¬ 
able abundances have been reported in the literature for most of 
the sources. In total we consider nine molecular species: SiO, 
H 2 S, c-HCCCH, HC 3 N, NH 2 CN, CH 3 CN, CH 3 OH, H 2 CO, and 
CH 3 CCH. Their abundances, together with the relevant refer¬ 
ences, are reported in Table lA.Tl For NGC 4418 the abundances 
were calculated assuming a molecular hydrogen colum n of 0.7- 
2.9x10^"* cm“^, as derived by ICostagliola et all (1201 3l) from CO 
2-1 observations. 

Our analysis includes well known Galactic objects, for which 
extensive literature is available, and six extragalactic sources 
for which mm/submm spectral scans have been performed. The 
Galactic objects span a wide range of physical conditions, in¬ 
cluding a hot-core (Sgr B2(N)), a typical nuclear bulge cloud 
(Sgr B2(OH)), the Orion Bar photo-dis sociation region, and the 
dark cloud TMC-1 (iMartm et al.ll200d and references therein). 
The extragalactic sources also s pan a wide range of galaxy types, 
including the ULIRG Arp 220 (iMartm et alJl201 ih, the evolved 
starb urst M 82, the young starbu rst NGC 253 (IMartm et al.l 
1200^ . and the Seyfert-2 NGC 1068 (lAladro et al.l2()13l). We also 
includ e the abundances derived in absorption bv Muller et"^ 
(I 2 OI ll) in the diffuse ISM of a z=0.89 galaxy located in the fore¬ 
ground of the blazar PKS 1830-211. 

The abundances and kinetic temperatures derived for each 
source are shown in Fig. | 6 ] A significant scatter in molecular 
properties is evident, even when considering the large uncertain¬ 
ties in the fit (color bar in the graph). Differences as high as four 
orders of magnitude are found, much higher than the typical dif¬ 
ferences in line ratios found for standard dense gas tracers, such 
as HCN, HNC and HCO*^, by extragalactic su r veys (factors of a 
few, e.g., iKrips et ^l20()8tlBaan et al]l2008t ICostagliola et al.l 
I 2 OIII) . 

It is evident that molecular abundances in NGC 4418 are 
very similar to those found in the much more massive and lumi¬ 
nous ULIRG Arp 220. The kinetic temperatures of most molec¬ 
ular species are of the order of 100 K, with complex molecules 
as CH 3 CN and CH 3 CCH having T^n >300 K. These tempera¬ 
tures are of the same order of magnitude of the ones found in hot 
cores and PDRs in the Galaxy, but in general higher than those 
found for the molecular gas of starburst galaxies. 


4.5.1. Principal component analysis 


In order to obtain a statistical description of the molecular prop¬ 
erties of th e sample, we performed a principal component anal¬ 
ysis (e.g., |Jolliffelll986l) . This is a standard procedure used to 
reduce the dimensionality of a problem and to hnd trends in 
multivariate data, and has been applied to extragalactic molec¬ 


ular data by various authors (e. g., Meier & Turneil I2005I: 


ICostagliola k al.ir201 it iMeier & Turneill2012l) . The uncertainty 
in the column density (N) and kinetic temperature (T) ht was 
taken into account by producing 1000 virtual samples for each 
source, with N and T values uniformly distributed inside their 


Article number, page 12 ofna 




































































































F. Costagliola et al.: Exploring the molecular chemistry and excitation in obscured luminous infrared galaxies 


Molecular abundances 























Fig. 6. Molecular abundances and kinetic te mperatures for Gala ctic and extragalactic objects. Data for Sgr B2(N), Sgr B2(OH), the Orion Bar, 
TMC-1, M 82, and NGC 253 were t aken from|Ma rtm_eta ] (|20()^, and references therein. Data f or Arp 220, NGC 1068, and PKS 1830-211 come 
from spectral scans, respectively bv iMartm et all (1201 ih . Aladro et al.) (l2013h . and iMuller et all (1201 ih . The color bars represent l-cr confidence 
intervals. In case no confidence interval was reported in the literature, a 20 % uncertainty was assumed. Filled bars represent upper limits. 
Abundances are reported in Table lA.7] 
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Fig. 7. Principal component analysis of molecular abundances. Left: Projection of the molecular abundances of Fig.[5]onto the first three principal 
components (PC). For each source, the l-tr confidence interval is shown as a shaded area. The color legend is the same as in Fig. [6] Right: Relative 
contribution of molecular abundances and kinetic temperatures to the first three principal components. For each PC, the upper and lower panels 
show contributions from abundance and temperature, respectively. The variables contributing for up to 75 % in total to each PC are highlighted in 
black. The fraction of the total variance (cr^) represented by each PC is shown in parentheses. 


confidence intervals. Assuming the abundance and kinetic tem¬ 
perature of all the analyzed molecular species as starting co¬ 
ordinates, an orthogonal set of principal com ponents (PC) was 
derive d following the procedure described in ICostagliola et al.l 
dMOl). As a result, we obtained the contribution of each observ¬ 
able to the nine PCs, and the projection of the sources on the PC 
basis. The uncertainties in the N and T fit result in a spread in 


the projected values, which is represented by the shaded areas in 

Fig.miefi). 

The first three PCs contain 80% of the variance in the sample, 
the contribution of each of the molecular properties being sum¬ 
marized in Fig. Ulfighf). The first principal component (PCI) 
is dominated by temperature and abundance of dense gas trac¬ 
ers, efficiently separating the cold gas of TMC-1 and the dense, 
warm environment of Arp 220 (see Fig.|7] upper left). The sec- 


Article number, page 13 ofM 

























































































































































































































A&A proofs: manuscript no. n4418_scan_v4 


ond component (PC2) is dominated by SiO, H 2 S, and c-HCCCH 
abundances, which appear to be anti-correlated with CH 3 CN and 
CH 3 OH abundances. Along this component, NGC 4418 has the 
highest projected value, while TMC-1 and the Orion Bar lie at 
the opposite end of the scale. The third component (PC3) is char¬ 
acterized by high HC 3 N and low CH 3 OH abundances, and high 
kinetic temperatures for all the considered tracers. This compo¬ 
nent efficiently separates NGC 4418 and Arp 220 from the rest 
of the Galactic and extragalactic objects. 

4.6. A new CON chemistry? 

Our PC analysis suggests that NGC 4418 and Arp 220 have 
molecular abundances and excitation which set them apart 
from all well-studied Galactic and extragalactic environments. 
Both objects are characterized by high kinetic temperatures 
(Tk >100 K) and high abundances of HC 3 N, SiO, H 2 S, and c- 
HCCCH. Compared to Galactic hot cores, the two sources show 
CH 3 OH abundances lower by more than one order of magnitude. 


Both galaxies have very compact (<100 pc) and warm 
(>100 K) molecular cores, and have been suggested to har¬ 
bor an obsc ured AGN and/or an extreme compac t star- 


burst (e.g.. 

Snoon et al. 2001bl 2004t Imanishi et al.l 

2004; 

Iwasawa et al. 

120051 

Downes & Eckartj 2007t Aalto et al. 

2009; 

Bateiat et al.l 

2011 ; 

Costasliola et al.l 120131; ISakamoto et al. 

20131 Barcos-Munoz et al. 2015). It is thus possible that the 


combination of strong radiation fields, extreme obscuration, and 
compactness of the molecular core are producing a peculiar 
chemistry which has not yet being investigated in detail. In¬ 
deed, current chemical models of PDRs, XDRs, and hot-cores 
struggle to reproduce the molecular abundanc es observed in the 
compact cores of LIRGs and ULIRGs (e.g., iKrins et ^120081 : 
ICostagliola et al.l [201 lUViti et alJl2014h . 

HCiN and c-HCCCH: Hot gas chemistry? In general, high 
HC 3 N abundances are associated with warm gas, shielded from 
UV radiation, typical of Galactic hot cores. In standard models, 
the molecule forms in the gas phase from acetylene and it is eas¬ 
ily de stroyed by ions and UV photons (C 2 H 2 - 1 -CN—»HC 3 N-i-H, 
e.g., lTurneretal.lh998h . This is in agreement w ith obser¬ 
vations of IC 342, where iMeier & Turneil (l2005h find that 
HC 3 N anti-correlates with PDRs. Howeve r, the molecule shows 
bright emission in dus t-obscured LIRGs dLindberg et alJl20Tll : 
ICostagliola et al.ll201 ih . which are thought to be powered either 
by a compact starburst or an AGN. Further more, recent obser- 
vations in the S eyfert galaxie s NGC 1097 dMartm et al.ll2015h 
and NGC 1068 dTakano et al.ll2014ll reveal bright HC 3 N emis¬ 
sion concentrated in the central 100 pc around the AGN (circum- 
nuclear disk, CND) and not in the more extended starburst ring. 
How the molecule can survive in such energetic environments is 
somewhat puzzling. A p ossible explan a tion m ay c ome from re¬ 
cent c hemical models bv iHarada et al.l d2010li and iHarada et al.l 
d2013h . which show that an enhancement of HC 3 N is possible 
in an AGN torus, thanks to hot (>300 K) gas phase chemistry. 
Abundances up to 10 ® can be reached in the plane of the disk, 
where the molecule is shielded from the ionizing radiation by 
large gas columns (i.e., N>10^"^ cm“^). Such a hot gas phase 
could also explain the high c- HCCCH abundance observed in 
NGC 4418 dHarada et al.ll2010h . 

CH^OH, SiO, and H 2 S: Shocks and XDRs ? Methanol forms 
on dust grains and is easily destr oyed by UV for A v<5 mag 
and by ions such as C^, He^ (e.g., iTurner et alJ[l998h . It is ob- 
seryed in the gas phase after eyaporation from dust mantles, and 
it is usually associated with young star formation (e.g.. Galactic 
hot-cores). The hot-core Sgr B2(N) is a clear example, haying a 


CH 3 OH abundance which is more than ten times higher than that 
found in all other sources in our sample (see Fig.|^. If the central 
100 pc of NGC 4 418 were dominated b y a young compact star- 
burst (<10 Myr, ISakamoto et al.ll^013h . we should obserye an 
enhancement of CH 3 OH when compared to more eyolyed star- 
bursts. The obseryed low abundance seems to be at odds with the 
young starburst scenario. 

Howeyer, the _ iKarby starburst galaxy Maffei 2, 

iMeier & Turneil d2012h find that CH 3 OH strongly corre¬ 
lates with SiO emission, and anti-correlates with sites of 
massiye star formation. The SiO molecule is liberated in the gas 
phase by the sputtering of dust particles due to strong shocks. 
The authors thus suggest that the two molecules may be tracing 
shocks rather than star formation. 

The presence of strong shocks in the nuclei of NGC 4418 and 
Arp 220 is supported by the high H 2 S abund ance. In their spec¬ 
tral scan with the SMA, iMartfn et alJ (1201 Ih find it to be higher 
in Arp 220 than in the starburtst NGC 253 by about two orders 
of magnitude. They interpret this oyerabundance to be the re¬ 
sult of grain disruption by shocks, and subsequent injection into 
the gas phase. This scenario is supported by the detection of a 
SiO ou tflow in the western nucleus of Arp 220 by lTunnard et al.l 
(I2015h . 

Since liberation of methanol from the grain mantles does not 
require sputtering, the relatiye abundances of SiO and CH 3 OH 
obseryed in NGC 4418 and Arp 220 may depend on the strength 
of the shocks. The two molecules haye yery different disso¬ 
ciation rates when exposed to ionizing radiation, the rate of 
CH 3 OH being more than te n times higher than for SiO (e.g., 
ISternberg & Dalgarnol Il995h . Therefore, an alternatiye expla¬ 
nation of the obseryed abundances may be the dissociation of 
CH 3 OH by X- or cosmic rays. 

In this scenario, HC 3 N should be also dissociated by the hard 
radiation and one would not expect such a high HC 3 N/CH 3 OH 
abundance ratio. Howeyer, methanol forms mainly on dust 
grains and has no efficient gas-phase formation route, while 
HC 3 N can be formed efficiently by hot gas phase reactions. The 
obseryed ratio may thus be a result of the different gas formation 
efficiencies when the dust grains haye been destroyed by shocks. 

In summary, the common chemistry of NGC 4418 and 
Arp 220 resembles that obseryed in Galactic hot-cores, but with 
a much lower methanol abundance and higher SiO. Young, dust- 
embedded star formation alone cannot explain the observed 
abundances and some kind of additional ionizing radiation is 
needed, together with mechanical feedback either from star for¬ 
mation or from an AGN. Dissociation of CH 3 OH may be due 
either to X-rays from an AGN or cosmic rays from super- 
noya remnants. If a starburst is powering the IR emission of 
NGC 4418, it must be in a young pre-supernoya stage (< 10 Myr 
ISakamoto et al.llT013h . and it is unclear whether the low super- 
noya rate could sustain the cosmic ray flux needed to dissoci¬ 
ate methanol. As an alternatiye, our obseryations are compatible 
with a composite starburst/AGN system in the inner 20 pc of the 
galaxy. 

5. Summary and Conclusions 

We report the results of a 70.7 GHz-wide spectral scan of the 
LIRG NGC 4418 with ALMA Cycle 0 obseryations in bands 3, 
6 and 7. Our spectral scan confirms that the chemical complexity 
in the nucleus of NGC 4418 is one of the highest eyer obseryed 
outside our Galaxy. We identify 317 emission lines aboye the 
three sigma leyel from a total of 45 molecular species, including 
15 isotopic substitutions and six yibrationally excited yariants. 
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We find that the molecular emission accounts for 15% of the 
total flux in band 3, and for 27% in bands 6 and 7. Our obser¬ 
vations clearly show that line contamination can have a serious 
impact on studies continuum emission in (U)LlRGs. This is es¬ 
pecially true for galaxies with broad lines, where the contami¬ 
nation from blended molecular emission could amount to up to 
30% of the flux. 

A combined LTE/NLTE fit of the spectrum reveals a multi¬ 
phase ISM, with temperatures ranging between 20 and more 
than 500 K, and densities between 10"^ and 10^ cm“^. The spec¬ 
trum of NCG 4418 is dominated by bright vibrationally excited 
HC 3 N, HCN , and HNC, with vibrational temperatures exceed¬ 
ing 350 K. We interpret this emission as the signature of an 
optically-thick, compact IR source (<5 pc of effective area) with 
brightness temperature >350 K in the core of the galaxy. We 
suggest that the presence of bright emission from vibrationally 
excited states of HC 3 N, HCN, and HNC may be a good tracer of 
dust-embedded compact sources in extragalactic objects. 

We compare the molecular abundances derived by our fit 
with those found in other Galactic and extragalactic environ¬ 
ments by means of a principal component analysis. We find 
that NGC 4418 and the ULIRG Arp 220 have similar molec¬ 
ular abundances and excitation, which set them apart from the 
other sources in the sample. We interpret this as the signature 
of a common, peculiar chemistry, which may be typical of the 
compact obscured nuclei of (U)LIRGs. 

While the line ratios of standard dense gas tracers such as 
HCN, HNC, and HCO^ can be reproduced by PDR models, 
other more complex molecules such as HC 3 N and c-HCCCH 
require either a hot-core chemistry or hot gas-phase reactions 
to be explained. The relative abundances of SiO and CH 3 OH 
seem to require a mix of strong shocks and dissociation by X- 
or cosmic-rays. The inferred chemistry, together with the strong 
vibrationally excited emission and compactness of the IR core, 
is consistent with a composite starburst/AGN system. 

Compared to Arp 220, NGC 4418 has a similarly bright and 
complex molecular spectrum, but much narrower line widths, 
which greatly reduce blending and simplify the identification of 
spectral features. The galaxy is thus the ideal target for future 
studies of molecular chemistry and excitation in IR-bright, ob¬ 
scured sources. 
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Appendix A: Fit detaiis for individuai moiecuies 

^^CO, C^^O: Carbon monoxide was detected in the J=l-0 tran¬ 
sitions of the two isotopic variants '^CO and C^^O. The latter 
is only a tentative detection, the peak flux being only at the 2- 
sigma level. Since we only detected one line per species, we 
could not fit an excitation temperature to the emission. By fixing 
the excitation temperature to 70 K, which is the lower limit to the 
CO_J=2-l brightness temperature detected bv ICostagliola et ^ 
( l2013h . a column density of the order of 10'* and 10*^ cm“^ is 
found for '*CO and C'*0, respectively. 

CS, '*CS, C**S, C*4S: We detect the J=2-l and J=6-5 ro¬ 
tational transitions of carbon monosulflde and its **C, **S, and 

isotopic variants. For CS, both the population diagram anal¬ 
ysis and the LTE fit And an excitation temperature of ~20 K, 
and a column density of 6+0.6x10*® cm“^. In order to calculate 
isotopic ratios, we fix the excitation temperature of the isotopic 
variants at 20 K. This implies the emission to be co-spatial and 
the excitation to be similar for all isotopic variants. Under this 
assumption, the LTE fit gives column densities of the order of 
10*® cm^^for '*CS and C**S, and of the order of lO'® cm^^for 
C*^S. An NLTE fit of CS emission was also performed. We And 
that the emission is best fit by a model with hydrogen density 
n(H2)=10® cm *, kinetic temperature 7'kin=50 K, and column 
density A=4xl0'^ cm“^. The NLTE values for Tkm and N are 
significantly higher than the LTE results. This may be due to 
the high critical density of the J=6-5 line, which would require 
n(H 2 )> 10 ^ cm~* to be efficiently excited by collisions. 

CN: We detect the CN multiplet at 3 mm, including nine 
merged lines spanning from 113.12 to 113.52 GHz in rest fre¬ 
quency. The lines have very similar upper-state energies (~5 K) 
and do not allow for a reliable fit. In order to obtain a lower 
limit to the molecule’s column density, we performed an LTE fit 
with a Axed temperature of 7’ex=70 K. The resulting column is 
N ~5xl0'^ cm“^. As discussed for '*CO and C'*0, the Tex value 
was chosen as the lower limit to the observed CO brightness 
temperature. 

NS: We detect the strongest lines of the 1=11/2-9/2 multi¬ 
plet of nitric sulflde at rest frequencies 253.57-255.6 GHz. Even 
if so close in frequency, the detected lines span a wide inter¬ 
val in upper state energies (£’u=40-360 K). We derive an ex¬ 
citation temperature of 7’ex=350±150 and a column density of 
A=8±4xl0*®cm“^. The large uncertainties on the derived quan¬ 
tities are due to the line blending in band 6. 

SO, *"*50: Emission from sulfur monoxide is heavily 
blended, the only sure detections being the lines 32-2i (99.3 
GHz), 23 -I 2 (109.25 GHz), and 65-55 (258.255). Because of the 
strong blending, the surface of our LTE At is not well de- 
flned, giving only a lower limit for the excitation temperature of 
7’ex>20 K. The fltted value for the column density is relatively 
well deflned, A=4±3xl0*® cm“^. We also performed a NLTE 
At of the isolated lines, resulting in hydrogen densities between 
10® and 10^ cm *, kinetic temperatures of 50-200 K, and col¬ 
umn densities from 2x10*® to 2x10*® cm“^. The isotopic vari¬ 
ant ®"*SO is only tentatively detected as two highly blended lines 
at 251.4 and 255.1 GHz. The LTE At with a Axed excitation tem¬ 
perature of 20 K results in a poorly constrained column density 
of 10*®-10*® cm-2. 

SiO, ^^SiO, ***810 : Strong silicon monoxide emission is 
detected in band 3 (J=2-l) and band 6 (J=5-4). The two rota¬ 
tional lines are free of blending and well identifled. Our LTE At 
gives an excitation temperature of 10+5 K and a column den¬ 
sity of 3+1x10*® cm“^, comparable with what found with the 
population diagram analysis. The NLTE At results in an hydro¬ 


gen density lower than 10® cm“®, a kinetic temperature greater 
than 20 K, and a column density greater than 10*® cm“^. We 
also detect the isotopic variants ^®SiO and ^'’SiO. The emission 
lines from these two species are strongly blended, but still clearly 
identifled. In our LTE At we Ax their excitation temperature to 
that of SiO, i.e. Tgx^lO K. The resulting column density estima¬ 
tion is of the order of 10*® cm“^, with large uncertainties (>50%) 
due to the line blending. 

HCN, H*®CN: Hydrogen cyanide is detected in its J=l-0 
and J=3-2 transitions. The J=3-2 line is the strongest molecu¬ 
lar line detected in our scan, with an integrated flux density of 
50+5 Jy km s *, and a peak brightness temperature of 35 K. Both 
emission lines have FWHM >140 kms *, broader than than the 
average 120 kms 'assumed in our At. The derived properties 
are thus representative only of the line core. The LTE At gives 
7’ex=7+l K and A=2+lxl0*® cm“^, while the NLTE analysis 
results in a hydrogen density lower than 10® cm“®, rkin> 100, 
and N >10*® cm“^. The large discrepancy in excitation temper¬ 
ature between the two fltting methods is most likely due to a 
limitation of the LTE At, which gives unreliable results for high 
opacities. An excitation temperature of 8 K seems too low when 
co mpared with the CO br ight ness temperature of >70 K reported 
bv ISakamoto et alJ (1201 3h and ICostagliola et aT] (1201 3h . and may 
be an indication of sub-thermal excitation of HCN. The bright¬ 
ness temperature for a source size of 0".4 is ~35 K and ~65 K 
for the two lines respectively (see Table |3. This sets a lower 
limit for the excitation temperature of the molecule, and thus to 
the kinetic temperature of the gas, of at least 65 K. Eor these 
reasons, we will assume the NLTE At to best represent the HCN 
emission in NGC 4418. 

The H*®CN J=l-0 line was also detected. A column density 
estimate of 2+1x10*® cm^^was obtained by Axing the excitation 
temperature at 7 K. If we assume an extr agalactic *^C/*®C ratio 
between 40 and 100 (iHenkel et al.ll2014h . this correspond to a 
column densities of HCN of (0.8-2)xl0*^ cm“^. 

HNC, HN*®C: We detect the J=l-0 and J=3-2 rotational 
transitions of hydrogen isocyanide. Both the population diagram 
analysis and the LTE At of the total spectrum result in an excita¬ 
tion temperature of 8-9 K, and column densities in the range 3- 
8x10*® cm“^. An NLTE At of the emission was also performed, 
resulting in an hydrogen density <10® cm ®, a kinetic tempera¬ 
ture greater than 100 K, and a column density >10*® cm-2. The 
results of the LTE and NLTE fltting of HCN and HNC are very 
similar, including the discrepancy of the fltted excitation tem¬ 
peratures for the two methods. The LTE At results in very low 
Tex, lower than the brightness temperature of the two detected 
lines of HNC (7’b=45 K and 30 K for the J=l-0 and J=3-2 tran¬ 
sitions, respectively). The NLTE At results in Tex-TO K, which 
is consistent with the brightness temperature of the lines for the 
assumed source size of 0".4. The J=l-0 transition of HN*®C is 
also detected. The LTE At with a Axed excitation temperature of 
8 K results in a column density of 1.5+0.5x10*® cm“^. 

HCN,v2=l, HNC,v2=l: Vibrationally excited HCN and 
HNC are tentatively detected in our scan. Eor each molecule 
we only detect the J=3-2,v2=lf transitions, therefore we can¬ 
not produce a reliable At to the rotational temperature. However, 
we can derive the vibrational temperature by comparing the in¬ 
tensity of the J=3-2 transitions of both the v2=l and v=0 states. 
We find that the relative population of the v2—l and v—0 states is 
well described by a vibrational temperature of ~350 Kfor HCN 
and 450 K for HNC. When assuming these temperatures to be 
also the rotational excitation temperatures of the v2=l emission, 
we obtain line intensities which are consistent with the observed 
spectrum. 
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HCO+, H'3C0+, The J=l-0 and J=3-2 transitions 

of HCO^ (formylium) are detected in band 3 and band 6 , respec¬ 
tively. The population diagram analysis and the LTE model fit 
result in an excitation temperature of 6-8 K and a column density 
of 1-6x10'^ cm“^. The RADEX models that best fits the emis¬ 
sion have n(H 2 )< 6 xlO^ cm“^, 7’kin>100 K, and N >10^® cm“^. 
Eor HCO^ we find the same discrepancy between the LTE and 
NLTE fitted excitation temperatures that we discuss in the HCN 
and HNC case. The NLTE fit provides more physically reason¬ 
able results, consistent with an observed brightness temperature 
of the lines of >20 K (see Table|2]i. 

We also detect the J=l-0 emission line of H'^CO^. 
By fixing the excitation temperature of the molecule at 
6 K, we obtain an LTE estimate of the column density of 
A(H^^CO^)=7+4xlO'"^ cm“^. If we assume an ext ragalactic 
ratio between 40 and 100 (iHenkel et al.ll2014l) . this cor¬ 
responds to a column density of HCO^ of ( 0 . 8 - 2 )xl 0 ^^ cm“^. 

The J=3-2 transition is tentatively detected in band 

6 , with heavy blending. The LTE fit is highly degenerate, and 
only results in an upper limit of 20 K for the excitation tempera¬ 
ture and a lower limit to the column density of lO'^cm^^. 

H 2 S: At least six transitions of hydrogen sulfide are included 
in the scan, but only one transition ( 2 ( 2 ,o)- 2 (i,i) at 216.71 GHz) 
is free of blending and clearly identified. By fixing the excita¬ 
tion temperature of H 2 S at 70 K, the LTE fit results in a column 
density of 1+0.5x10^^ cm“^. 

CCH: We detect the N=l-0 multiplet of ethynyl at rest fre¬ 
quencies between 87.28 and 87.45 GHz (sky frequency ~86.76 
GHz). The six lines have very similar upper state energies (~4 K) 
and thus do not allow for an accurate fit of the excitation tem¬ 
perature. By fixing T^^-IQ K we obtain an LTE estimate of the 
column density of 2 + 1 x 10 *^ cm“^. 

HCS^: The J=2-l and J=6-5 rotational transitions of HCS^ 
(thioformylium) are detected in band 3 and band 6 , respectively. 
The 1 mm line is heavily blended with emission from vibra- 
tionally excited HC 3 N, while the 3 mm line is free of blend¬ 
ing. The emission can be reproduced by our LTE model with 
rex=20+10 K and A=8+4xl0'^. 

CCS: We detect six emission lines of thioxoethenylidene, 
with upper energy levels ranging from 20 to 110 K. The emission 
is free of blending and is well fit by our LTE model. The best fit 
values are 7’ex=20+10 K and 1.5+0.5x10'® cm“^. 

N 2 H^ : We detect the J=l-0 triplet at 93.17 GHz and the 
J=3-2 line of N 2 H^ (diazenylium). The triplet emission is totally 
blended, but free of contamination from other species. The LTE 
and NLTE fit procedures give very similar results (see Table |4li 
with 7’ex~30 K and N ~5xl0'® cm“^. The best fit hydrogen den¬ 
sity in our NLTE analysis is lO'^ cm“^, which is higher than the 
critical density of the detected transitions (mc ~10®-10® cm“^), 
and the Tkin and Tex coincide for both lines. We conclude that the 
RADEX best fit represents an LTE excitation of the molecule. 

H 2 CO: We detect eight transitions of formaldehyde, four 
from para-H 2 CO and four from ortho-H 2 CO. The LTE fit for 
both species results in an excitation temperature of 350+100 K 
and column a column density of 5-8x10'® cm“^. The NLTE anal¬ 
ysis results in an hydrogen density of 10 ® cm“®, a kinetic tem¬ 
perature greater than 300 K and a column density of 7x10'®- 
10 '® cm“^. In our NLTE fit we only obtain a lower limit to the 
kinetic temperature because the LAMBDA database does not 
contain collision coefficients for temperatures greater than 300 
K. 

c-HCCCH: We detect 15 emission lines of cyclopropenyli- 
dene at the 3-cr level across ALMA bands 3 and 6 . The de¬ 
tected lines have upper state energies ranging from 40 to 500 K. 


The LTE model well fits the emission, resulting in an exci¬ 
tation temperature of 180+100 K and a column density of 
1+0.5x10'^ cm“^. The large uncertainties are mostly due to the 
blending in band 6 . 

H 2 CS, H}^CS: Three isolated emission lines of thio- 
formaldehyde are detected at sky frequencies 103.89, 268.66, 
and 276.95 GHz. Three more lines from Hj^CS are tentatively 
detected (~ 2cr) in band 3. The detected transitions of H 2 CS 
have upper state energies ranging from 20 to 70 K. The LTE fit 
results in an excitation temperature of 35+10 K and a column 
density of 1+0.5x10'® cm“^. By fixing the excitation tempera¬ 
ture at 35 K we derive an upper limit for the column density of 
H^®CS of 5+2x10'® cm-2. 

HC 3 N, H'®CCCN, HCC'®CN: We detect bright 
cyanoacetylene emission in all three ALMA bands. Nine 
isolated lines are detected, with upper state angular momentum 
ranging from 10 to 32, and upper state energies from 24 to 
230 K (see Table |2]). The population diagram analysis and the 
LTE fit of the whole spectrum give an excitation temperature 
of ~70 k and column densities in the range 2-9x10'® cm“^(see 
Tables |2]|4l). The H'®CCCN and HCC'®CN isotopomers are also 
detected in all three bands, for a total of nine isolated emission 
lines with upper state energies ranging from 20 to 240 K. 
We performed an LTE fit of both species, assuming the same 
excitation temperature of 70 K as HC 3 N. The fitted column 
densities are very similar, in the range of 1 - 2 x 10 '® cm“^, with 
uncertainties of the order of 50%. 

Vibrationally excited HC 3 N: Emission from vibrationally 
excited HC 3 N dominates the mm-wave spectrum of NGC 4418. 
We detect rotational transitions from the four vibrationally ex¬ 
cited states HC 3 N,v 6 = 1 , HC 3 N,v 7 = 1 , HC 3 N,v 6 = 1 ,v 7 = 1 , and 
HC 3 N,v 7 = 2 , for a total of 25 transitions. The upper state en¬ 
ergies of such vibrationally excited lines are very high, ranging 
from 300 to 900 K. The emission from all four species is well 
fitted by our LTE model and by the rotational diagram analysis. 
We find similar excitation temperatures, in the range 70-100 K, 
and column densities ranging from 10'®-10'® cm“^. Erom the 
relative population of the vibrational levels it is possible to de¬ 
rive the vibrational temperature of the molecule via a population 
diagram analysis. The results are reported in Eig.|4l We are able 
to derive the vibrational temperature of three different angular 
momentum states (J= 10,11,12) for which emission of more than 
two vibrational states was detected. The three fits give remark¬ 
ably similar results, with vibrational temperatures ranging from 
330 to 400 K. 

CH 2 NH: We identify three spectral features as emission 
from methanimine at rest frequencies 105.79, 257.11, and 
266.27 GHz, with upper state energies ranging from 30 to 60 K. 
The LTE fit of the spectrum results in an excitation temperature 
of 45+30 K and a column density of 1.5+0.5x10'^ cm“^. 

NH 2 CN: Emission from cyanamide is tentatively detected 
in NGC 4418. Only one isolated line at sky frequency 98.62 
GHz is unambiguously identified as the NH 2 CN J=5-4 transi¬ 
tion, while most emission across the observed ALMA bands is 
heavily blended. Our LTE fit results in an excitation temperature 
of 30+20 K and a column density of 5+2x10'® cm“^. However, 
the derived parameters are highly uncertain because of the blend¬ 
ing of the lines at 1 mm. 

CH 3 CN: Bright emission from methyl cyanide is detected as 
four line forests at sky frequencies 91.3-91.35, 109.56-109.62, 
254.87-255.78,292.11-292.28 GHz. The detected lines have up¬ 
per state energies ranging from 13 to 500 K, with bright emission 
from lines with Eup >200 K. The LTE fit results in a very high 
excitation temperature of 1000 K, but does not reproduce the 
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intensity of the lines at 3 mm. The NLTE procedure provides a 
better fit of the emission across all three ALMA bands, resulting 
in a low hydrogen density of 10 "^ cm“^, and a high kinetic tem¬ 
perature >500 K, required in order to fit the high energy lines. 

The critical density of the observed transitions varies between 
10^ cm“^ for the lines at 3 mm, and 10^ cm“^ for the lines at 
1 mm. It appears that the only way to fit both the 3 mm and the 
high-energy 1 mm lines is to lower the hydrogen density until the 
3 mm lines are not efficiently excited by collisions («(H 2 )< n^, 
i.e. they are not at LTE). This could imply that CH 3 CN emission 
is originating from a gas phase that is significantly more diffuse 
than for the other molecular tracers, or that the low- and high- 
excitation lines are emerging from different environments. 

CH^CCH: Methyl acetylene is tentatively detected in 
NGC 4418. The emission is heavily bended and no isolated line 
can be identified. Emission from the molecule is required in or¬ 
der to fit the line forests at sky frequencies 254-254.6, 271.2- 
271.6, and 288-288.6 GHz. The LTE fit results in an excitation 
temperature of 350+100 K, and a column density of 1 +0.5x10^®. 

We advise extreme caution when interpreting these results, since 
the strong blending makes the fit highly degenerate. 

HC 5 N: Emission from cyanobutadiyne is required in order to 
explain 11 spectral features across all three ALMA bands, which 
could not be associated to any other molecule. Nine of the de¬ 
tected lines are isolated lines free of blending with other species. 

The upper state angular momentum of the detected transitions 
ranges from 32 to 42, corresponding to energies above ground 
of 65 to 85 K. The population analysis and the LTE fit of the 
whole spectrum result in an excitation temperature of 60-70 K, 
and column densities in the range 0.5-2x10'® cm“^. 

CH^OH: Methanol emission is needed in order to fit the 
highly blended region between sky frequencies 287.8 and 288.3 
GHz. No other transition from the molecule is detected above 
the 3-cr level across the whole observed frequency range. The 
NLTE fit results in a kinetic temperature of 100+60 K and a 
column density of 2 ±lxl 0 '® cm“^, for an hydrogen density of 
10^ cm“^. The fit is poorly constrained because of line confusion. 

Non-LTE excitation of methyl acetylene could give a significant 
contribution to the emission in this region, thus the derived col¬ 
umn density of CH 3 OH should be considered as an upper limit. 


Article number, page 18 ofl34l 


F. Costagliola et al.: Exploring the molecular chemistry and excitation in obscured luminous infrared galaxies 

Rest Frequency [GHz] 


84.57 84.87 85.18 85.48 85.79 86.09 86.40 86.70 



86.50 86.79 87.07 87.36 87.65 87.94 88.23 88.51 



88.31 88.62 88.92 89.23 89.53 89.84 90.14 90.45 



90.24 90.60 90.96 91.32 91.67 92.03 92.39 92.74 



Sky Frequency [GHz] 


Fig. A.l. The ALMA spectral scan of NGC 4418. The data are shown as a shaded histogram plot, the best fit for each molecule is shown as solid 
black lines. The total fit to the data is shown as a red solid line. The l-cr level is shown as a dashed green line. The residuals, together with the data 
rms, are shown below the spectrum as a green histogram plot with error bars. Emission lines above 3-(T are labeled. 
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Fig. A.2. Continues from Fig. IA.il 
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Fig. A.3. Continues from Fig. IA.il 
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Fig. A.4. Continues from Fig. IA.il 
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Fig. A.5. Continues from Fig. lA. 1 
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Fig. A.6. Continues from Fig. IA.il 
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Fig. A.7. Continues from Fig. IA.il 
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Fig. A.8. Continues from Fig. lA. 1 
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Fig. A.9. Continues from Fig. IA.il 
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Table A.l. Molecular lines detected above the 3-sigma level. A total of 317 lines were detected in our spectral scan. The peak flux densities 
reported in table were derived with the LTE/NLTE fit of the whole spectrum. The signal to noise ratio reported in column six is the ratio between 
the peak flux density and the average rms of the ALMA band containing the line. The notation for the transition quantum numbers is the one used 
by the Splatalogue database. Tentative detections are shown in italics. 
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Sky freq. [GHz] 
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Table A.2. Continues from Table lA. l1 
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Sky freq. [GHz] 

Molecule 

Transition 

Rest freq. [GHz] 

Peak [mJy] 
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250.082 

SO 

6;5-5;4 

251.825 

55 

67 

250.123 

CHiOH 

4;-3;0-4;2;0 

251.867 

3 

3 

250.147 

CHiOH 

5;3;0-5;-2;0 

251.891 

3 

4 

250.152 

CHiOH 

6;3;0-6;-2;0 

251.896 

3 

3 

250.157 

CH^OH 

4;3;0-4;-2;0 

251.900 

3 

3 

250.390 

c-HCCCH 

9(2,7)-9(l,8) 

252.136 

8 

9 

250.390 

c-HCCCH 

9(3,7)-9(2,8) 

252.136 

3 

3 

250.662 

c-HCCCH 

8(2,7)-8(l,8) 

252.410 

4 

5 

251.344 

c-HCCCH 

22(20,3)-22(19,4) 

253.096 

3 

4 

251.454 


6(6)-5(5) 

253.207 

14 

17 

251.815 

NS 

J=ll/2-9/2,Q=l/2,F=13/2-ll/2,l=e 

253.570 

15 

18 

251.815 

NS 

J=1 l/2-9/2,0=l/2,F=l l/2-9/2,l=e 

253.570 

12 

15 

251.817 

NS 

J=1 l/2-9/2,Q=l/2,F=9/2-7/2,l=e 

253.572 

10 

12 

251.887 

HCC'^CN 

J=28-27 

253.644 

8 

10 

252.210 

NS 

J=ll/2-9/2,Q=l/2,F=13/2-ll/2,l=f 

253.968 

15 

18 

252.212 

NS 

J=1 l/2-9/2,Q=l/2,F=l l/2-9/2,l=f 

253.971 

12 

15 

252.212 

NS 

J=1 l/2-9/2,Q=l/2,F=9/2-7/2,l=f 

253.971 

10 

12 

252.457 

3“SiO 

6-5 

254.217 

34 

41 

252.774 

c-HCCCH 

24(21,4)-24(20,5) 

254.536 

2 

3 

252.800 

c-HCCCH 

25(17,8)-25(16,9) 

254.563 

3 

4 

252.922 

CH 2 NH 

4(0,4)-3(0,3) 

254.685 

89 

109 

252.936 

HC 3 N 

J=28-27 

254.700 

118 

144 

253.222 

c-HCCCH 

5(3,3)-4(2,2) 

254.988 

4 

5 

253.351 

HC 3 N v 6 =l 

J=28-27,l=le 

255.117 

26 

32 

253.550 

HC 3 N v 6 =l 

J=28-27,l=lf 

255.317 

26 

32 

253.557 

HC 3 N v7=l 

J=28-27,l=le 

255.325 

79 

97 

253.711 

HC'^ 6 )+ 

3-2 

255.479 

22 

27 

253.827 

NS 

J=ll/2-9/2,0=3/2,F=13/2-ll/2,l=f 

255.597 

6 

7 

253.827 

NS 

J=ll/2-9/2,0=3/2,F=13/2-ll/2,l=e 

255.597 

6 

7 

253.831 

NS 

J=1 l/2-9/2,Q=3/2,F=l l/2-9/2,l=f 

255.600 

5 

6 

253.831 

NS 

J=1 l/2-9/2,Q=3/2,F=l l/2-9/2,l=e 

255.600 

5 

6 

253.833 

NS 

J=1 l/2-9/2,Q=3/2,F=9/2-7/2,l=f 

255.603 

4 

5 

253.833 

NS 

J=1 l/2-9/2,0=3/2,F=9/2-7/2,l=e 

255.603 

4 

5 

253.870 

H'^CCCN 

J=29-28 

255.640 

15 

19 

253.919 

HC 3 N v7=l 

J=28-27,l=lf 

255.689 

79 

97 

254.069 

CH 2 NH 

4(2,3)-3(2,2) 

255.840 

37 

46 

254.175 

HC 3 N v6=l,v7=l 

(28 0 0)-(7 0 0) 

255.947 

10 

12 

254.216 

HC 3 N v6=l,v7=l 

(28 0 l)-(7 0 1) 

255.988 

9 

12 

254.255 

HCS+ 

6-5 

256.027 

35 

43 

254.299 

HC 3 N v6=l,v7=l 

(28-2 2)-(7 2 2) 

256.072 

9 

11 


Table A.3. Continues from Table lA. f] 
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Sky freq. [GHz] 

Molecule 

Transition 

Rest freq. [GHz] 

Peak [mJy] 

SNR 

254.338 

HC 3 N v6=l,v7=l 

(28 2 2)-(7-2 2) 

256.111 

9 

12 

254.387 

CH 3 CCH 

15(6)-14(6) 

256.161 

5 

6 

254.392 

CHjNH 

4(3,2)-3(3,l) 

256.165 

9 

11 

254.403 

CH 2 NH 

4(3,l)-3(3,0) 

256.177 

9 

11 

254.441 

CH 3 CCH 

15(5)-14(5) 

256.214 

3 

4 

254.484 

CH 3 CCH 

15(4)-14(4) 

256.258 

4 

5 

254.485 

HC 3 N v7=2 

J=28-27,l=0 

256.260 

32 

39 

254.518 

CH 3 CCH 

15(3)-14(3) 

256.293 

10 

12 

254.537 

HC 3 N v7=2 

J=28-27,l=2e 

256.311 

31 

37 

254.543 

CH 3 CCH 

15(2)-14(2) 

256.317 

6 

7 

254.557 

CH 3 CCH 

15(1)-14(1) 

256.332 

6 

7 

254.562 

CH 3 CCH 

15(0)-14(0) 

256.337 

6 

8 

254.591 

HC 3 N v7=2 

J=28-27,l=2f 

256.366 

31 

37 

255.034 

CH 3 CN 

14;12-13;12 

256.812 

4 

4 

255.099 


7(6)-6(5) 

256.878 

19 

23 

255.148 

CH 3 CN 

14;11-13;11 

256.926 

4 

5 

255.251 

CH 3 CN 

14;10-13;10 

257.031 

4 

5 

255.333 

CH 2 NH 

4(2,2)-3(2,l) 

257.113 

38 

46 

255.345 

CH 3 CN 

14;9-13;9 

257.125 

12 

15 

255.429 

CH 3 CN 

14;8-13;8 

257.209 

6 

7 

255.474 

25SiO 

6-5 

257.255 

26 

32 

255.503 

CH 3 CN 

14;7-13;7 

257.284 

7 

8 

255.567 

CH 3 CN 

14;6-13;6 

257.349 

17 

21 

255.621 

CH 3 CN 

14;5-13;5 

257.403 

8 

10 

255.666 

CH 3 CN 

14;4-13;4 

257.448 

9 

11 

255.700 

CH 3 CN 

14;3-13;3 

257.483 

22 

27 

255.725 

CH 3 CN 

14;2-13;2 

257.507 

10 

12 

255.740 

CH 3 CN 

14;1-13;1 

257.522 

10 

13 

255.744 

CH 3 CN 

14;0-13;0 

257.527 

12 

15 

256.381 

NH 2 CN 

13(1,13)-12(1,12), 

258.169 

15 

18 

256.467 

SO 

6;6-5;5 

258.255 

63 

77 

261.966 

HC 3 N 

J=29-28 

263.792 

112 

137 

262.158 

CH 2 NH 

3(2,1)-4(1,4) 

263.986 

8 

10 

262.162 

H213CS 

8(0,8)-7(0,7) 

263.989 

5 

7 

262.395 

HC 3 N v 6 =l 

J=29-28,l=le 

264.224 

27 

33 

262.601 

HC 3 N v 6 =l 

J=29-28,l=lf 

264.432 

27 

33 

262.609 

HC 3 N v7=l 

J=29-28,l=le 

264.440 

79 

97 

262.621 

Hi^CCCN 

J=30-29 

264.451 

15 

18 

262.636 

c-HCCCH 

21(20,1)-21(19,2) 

264.467 

3 

3 

262.984 

HC 3 N v7=l 

J=29-28,l=lf 

264.817 

80 

97 

263.244 

HC 3 N v 6 =l, v7=l 

(29 0 0)-(8 0 0) 

265.079 

10 

12 

263.289 

HC 3 N v 6 =l, v7=l 

(29 0 l )-(8 0 1) 

265.125 

10 

12 

263.380 

HC 3 N v 6 =l, v7=l 

(29 2 2)-(8-2 2) 

265.216 

9 

12 

263.424 

HC 3 N v 6 =l, v7=l 

(29-2 2)-(8 2 2) 

265.260 

9 

12 

263.566 

HC 3 N v7=2 

J=29-28,l=0 

265.404 

32 

39 

263.920 

c-HCCCH 

4(4,l)-3(3,0) 

265.759 

16 

19 

264.012 

HCN v2=l 

J=3-2,l=le 

265.853 

47 

58 

264.046 

HCN 

3-2 

265.886 

379 

465 

264.347 

c-HCCCH 

23(15,8)-23(14,9) 

266.190 

5 

6 

264.427 

CH 2 NH 

4(1,3)-3(1,2) 

266.270 

80 

98 

265.020 

c-HCCCH 

24(17,8)-24(16,9) 

266.867 

4 

5 

265.195 

c-HCCCH 

23(21,2)-23(20,3) 

267.044 

3 

3 

265.349 

HCN v2=l 

J=3-2,l=lf 

267.199 

48 

59 

265.705 

HCO+ 

3-2 

267.558 

213 

261 

266.157 

m^^cs 

8(1,7)-7(1,6) 

268.012 

11 

14 

268.648 

H 2 CS 

8(1,8)-7(1,7) 

270.521 

22 

27 


Table A.4. Continues from Table lA. l1 
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Sky freq. [GHz] 

Molecule 

Transition 

Rest freq. [GHz] 

Peak [mJy] 

SNR 

269.872 

HCC'^CN 

J=30-29 

271.753 

8 

9 

269.914 

c-HCCCH 

22(15,8)-22(14,9) 

271.795 

7 

8 

270.042 

HNC v2=l 

J=3-2,l=le 

271.924 

86 

53 

270.098 

HNC 

3-2 

271.981 

328 

201 

270.863 


6-5 

272.752 

46 

28 

270.996 

HC 3 N 

J=30-29 

272.885 

104 

64 

271.341 

CH 3 CCH 

16(6)-15(6) 

273.232 

6 

4 

271.344 

c-HCCCH 

21(13,8)-21(12,9) 

273.236 

8 

5 

271.371 

Hi^CCCN 

J=31-30 

273.263 

14 

9 

271.439 

HC 3 N v 6 =l 

J=30-29,l=le 

273.332 

27 

16 

271.444 

CH 3 CCH 

16(4)-15(4) 

273.337 

5 

3 

271.480 

CH 3 CCH 

16(3)-15(3) 

273.373 

12 

8 

271.506 

CH 3 CCH 

16(2)-15(2) 

273.399 

7 

4 

271.522 

CH 3 CCH 

16(1)-15(1) 

273.415 

8 

5 

271.527 

CH 3 CCH 

16(0)-15(0) 

273.420 

8 

5 

271.653 

HC 3 N v 6 =l 

J=30-29,l=lf 

273.546 

27 

16 

271.660 

HC 3 N v7=l 

J=30-29,l=le 

273.554 

79 

48 

275.535 

i^CS 

6-5 

277.455 

25 

15 

275.853 

c-HCCCH 

19(11,8)-19(10,9) 

277.776 

11 

7 

276.095 

NH 2 CN 

14(1,14)-13(1,13), 

278.020 

12 

8 

276.714 

CH 2 NH 

8(l,7)-8(0,8) 

278.643 

29 

18 

276.956 

H 2 CS 

8(1,7)-7(1,6) 

278.886 

23 

14 

277.577 

N 2 H+ 

3-2 

279.512 

144 

88 

277.604 

c-HCCCH 

18(11,8)-18(10,9) 

279.539 

12 

8 

277.882 

NH 2 CN 

14(0,14)-13(0,13), 

279.819 

7 

4 

278.864 

HCC^^CN 

J=31-30 

280.808 

7 

4 

278.922 

c-HCCCH 

17(9,8)-17(8,9) 

280.867 

14 

8 

279.578 

0 -H 2 CO 

4;1;4-3;1;3 

281.527 

100 

61 

279.758 

NH 2 CN 

14(1,13)-13(1,12), 

281.708 

12 

8 

280.025 

HC 3 N 

J=31-30 

281.977 

97 

59 

280.085 

c-HCCCH 

16(8,8)-16(7,9) 

282.037 

5 

3 

280.096 

c-HCCCH 

16(9,8)-16(8,9) 

282.049 

15 

9 

280.121 

Hi^CCCN 

J=32-31 

282.073 

13 

8 

280.426 

c-HCCCH 

4(4,0)-3(3,l) 

282.381 

5 

3 

280.483 

HC 3 N v 6 =l 

J=31-30,l=le 

282.439 

27 

16 

280.703 

HC 3 N v 6 =l 

J=31-30,l=lf 

282.660 

27 

16 

280.711 

HC 3 N v7=l 

1=31-30,l=le 

282.668 

77 

47 

281.053 

c-HCCCH 

15(7,8)-15(6,9) 

283.012 

16 

10 

281.056 

c-HCCCH 

15(8,8)-15(7,9) 

283.016 

5 

3 

281.112 

HC 3 N v7=l 

J=31-30,l=lf 

283.072 

78 

47 

281.378 

HC 3 N v 6 =l, v7=l 

(31 OO)-(OOO) 

283.340 

10 

6 

281.434 

HC 3 N v 6 =l, v7=l 

(31 0 l)-(0 0 1) 

283.396 

9 

6 

281.539 

HC 3 N v 6 =l, v7=l 

(31 2 2)-(0-2 2) 

283.502 

9 

6 

281.594 

HC 3 N v 6 =l, v7=l 

(31-2 2)-(0 2 2 ) 

283.557 

9 

6 

287.207 

c^^s 

6-5 

289.209 

56 

34 

287.855 

HCC^^CN 

J=32-31 

289.862 

7 

4 

288.102 

CH^OH 

6;0;0-5;0;0 

290.111 

13 

8 

288.293 

CH 3 CCH 

17(6)-16(6) 

290.303 

8 

5 

288.353 

CH 3 CCH 

17(5)-16(5) 

290.364 

5 

3 

288.403 

CH 3 CCH 

17(4)-16(4) 

290.413 

6 

4 

288.441 

CH 3 CCH 

17(3)-16(3) 

290.452 

15 

9 

288.469 

CH 3 CCH 

17(2)-16(2) 

290.480 

9 

5 

288.485 

CH 3 CCH 

17(1)-16(1) 

290.497 

9 

6 

288.491 

CH 3 CCH 

17(0)-16(0) 

290.502 

9 

6 

288.551 


6(7)-5(6) 

290.562 

19 

12 

288.611 

P-H 2 CO 

4;0;4-3;0;3 

290.623 

63 

39 


Table A.5. Continues from Table lA. f] 
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Sky freq. [GHz] 

Molecule 

Transition 

Rest freq. [GHz] 

Peak [mJy] 

SNR 

288.870 

H'^CCCN 

J=33-32 

290.884 

12 

7 

289.053 

HC 3 N 

J=32-31 

291.068 

88 

54 

289.221 

P-H 2 CO 

4;2;3-3;2;2 

291.238 

40 

24 

289.363 

0 -H 2 CO 

4;3;2-3;3;l 

291.380 

34 

21 

289.367 

0 -H 2 CO 

4;3;l-3;3;0 

291.384 

33 

20 

289.468 


6-5 

291.486 

20 

12 

289.527 

HC 3 N v 6 =l 

J=32-31,l=le 

291.545 

26 

16 

289.754 

HC 3 N v 6 =l 

J=32-31,l=lf 

291.774 

26 

16 

289.762 

HC 3 N v7=l 

J=32-31,l=le 

291.782 

75 

46 

289.927 

P-H 2 CO 

4;2;2-3;2;l 

291.948 

36 

22 

290.176 

HC 3 N v7=l 

J=32-31,l=lf 

292.199 

76 

46 

290.444 

HC 3 N v 6 =l, v7=l 

(32 0 0)-(l 0 0) 

292.469 

9 

6 

290.506 

HC 3 N v 6 =l, v7=l 

(32 0 1)-(1 0 1) 

292.531 

9 

6 

290.619 

HC 3 N v 6 =l, v7=l 

(32-2 2)-(l 2 2) 

292.645 

9 

6 

290.647 

CH^OH 

6;-l;0-5;-l;0 

292.673 

8 

5 

290.679 

HC 3 N v 6 =l, v7=l 

(32 2 2)-(l-2 2) 

292.705 

9 

6 

290.759 

H2'3cS 

9(1,9)-8(1,8) 

292.786 

11 

7 

290.804 

HC 3 N v7=2 

J=32-31,1=0 

292.832 

31 

19 

290.881 

HC 3 N v7=2 

J=32-31,l=2e 

292.909 

30 

19 

290.961 

HC 3 N v7=2 

J=32-31,l=2f 

292.990 

30 

19 

291.808 

CH 3 CN 

16;9-15;9 

293.843 

8 

5 

291.877 

CS 

6-5 

293.912 

251 

154 

292.062 

CH 3 CN 

16;6-15;6 

294.098 

12 

7 

292.124 

CH 3 CN 

16;5-15;5 

294.161 

6 

3 

292.175 

CH 3 CN 

16;4-15;4 

294.212 

6 

4 

292.214 

CH 3 CN 

16;3-15;3 

294.251 

15 

9 

292.242 

CH 3 CN 

16;2-15;2 

294.280 

7 

4 

292.259 

CH 3 CN 

16;1-15;1 

294.297 

7 

4 

292.265 

CH 3 CN 

16;0-15;0 

294.302 

8 

5 


Table A.6. Continues from Table lA. f] 


NGC4418 Arp220 
SiO -9.2...-8.2 -9.4...-8.9 
H 2 S -7. 8 ...-6.7 -7.0...-6.7 
c-HCCCH -7.8...-6.7 - 8 . 6 ...-8.2 
HC 3 N -8.5...-7.2 -8.0...-7.7 
NH 2 CN -9.0...-8.0 
CH 3 CN -8.4...-7.6 -8.5...-8.1 
CH 3 OH -8.5...-7.4 -1.1...-1.2 
H 2 CO - 8 .3...-7.5 -8.9...-8.5 
CH 3 CCH -7. 8 ...-6.7 -7.5...-7.1 


M82 

NGC253 

NGC 1068 

- 


- 10 . 2 .. 

.- 10.0 

-9.5.. 

.-9.3 

-9.6.. 

.-8.5 

-9.5.. 

.-9.1 

- 


-9.3.. 

.- 8.8 

-9.4.. 

.-9.1 

- 10 . 1 . 

..-9.8 

-9.7.. 

.-8.9 

-9.7.. 

.-9.3 

-8.9.. 

.- 8.1 

-9.9.. 

.-9.0 

- 10 . 0 . 

..-9.5 

<-9.7 

- 10 . 2 . 

..-9.6 

-9.7.. 

.-9.3 

-9.6.. 

.-9.5 

-9.6.. 

.- 8.2 

-7.9.. 

.-7.7 

-7.9.. 

.-7.7 

-8.9.. 

.-8.4 

-9.O.. 

.-8.7 

<-l 

M 

-7.8.. 

.-7.4 

- 8 . 2 .. 

.- 8.0 

<-8.5 


SgrB2(N) 

SgrB2(OH) 

Orion 

Bar 

TMC-l 

PKS 1830 

-10.9... 

-10.3 

- 


-10.4... 

.- 10.2 

<-l 

1.6 

-9.5.. 

.-9.3 

- 10 . 2 .. 

.-9.5 

- 


-8.3... 

.- 8.1 

<-S 

1.3 

- 


- 11 . 0 .. 

.-9.9 

-10.4.. 

..-9.5 

-9.8... 

-9.6 

-9.4.. 

.-9.2 

-8.7.., 

.-8.5 

-7.8... 

-7.0 

-9.5... 

-8.9 

- 


-8.3.. 

.- 8.1 

-9.4.. 

.-9.2 

- 11 . 1 .. 

.-9.6 

-10.4.. 

..-9.4 

- 


- 


- 


-7.0... 

- 6.2 

-9.9... 

-9.3 

<-10.3 

-9.1.. 

.-8.9 

-9.3.. 

.-9.1 

- 6 . 1 ... 

-5.4 

-8.3... 

.-7.2 

<- 10.1 

- 8 . 8 .. 

.- 8.6 

- 8 . 2 .., 

.- 8.0 

-9.6... 

-8.9 

-8.3... 

.-7.6 

-8.3... 

.- 8.1 

-7.8.. 

.-7.6 

- 8 . 0 .. 

.-7.8 

- 8 . 6 ... 

-7.9 

-9.1... 

.-8.4 

<-8 

;.5 

-8.3.. 

.- 8.1 

- 8 . 6 .. 

.-8.4 


Table A.7. Base-ten logarithm of molecular abundances in Galactic and extragalactic sources. For NGC 4418, the range of values represents the 
confidence interval of the m inimization in the best fit model. Data for Sgr B2(N), Sgr B2(OFt) , the Orion Bar, TMC-1, and NGC 253 were 
taken fromiMartm et al.l ( [^06h . Data for Arp 220, NGC 1068, and PKS 1830-211, were taken from iMartm et al.l ( 1201 Ih . lAladro et alj ( l2013h . and 
iMuller et alJ ll201 ill respectively. 
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HCN 1-0 


HCO‘^l-0 


HNC 1-0 


0.08 
g 0.06 

■t 0.04 

= 

0.02 
E 0 


A 

a 0.04 



0.06 

^ 0.04 
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Fig. A.IO. Comparison of our ALMA data with single-dish spectra of the brightest lines in the scan. The blue line shows the spectrum extracted 
from the ALMA visibilities by fitting a point source with the CASA routine uvmultifit. The black line shows single dish spectra observed by 
the authors with the IRAM 30m telescope. The error bars show the squared sum of the rms of the two datasets. The spectra were interpolated to a 
common resolution of 28 and 19.5 MHz at 3 and 1 mm, respectively. 
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Fig. A.ll. Comparison of our ALMA data with SMA observations bv ICostagliola et aP l l2013h . The blue line shows the spectrum extracted from 
the ALMA visibilities by fitting a point source with the CASA routine uvmultifit. The black line shows the flux density measured by the 
extended configuration of the SMA, with a beam size of 0".4 and a maximum recoverable scale of 2". The error bars show the squared sum of the 
rms of the two datasets. The spectra were interpolated to the resolution of the ALMA spectral scan. 


Table A.8. Integrated flux densities from pre-ALMA observations 


Molecule 


Integrated flux densities [Jy km/s] 



ALMA'‘> 

IRAM 30m'^> 

SMA''> 

HCN 1-0 

11 .2+1 

13+2 

_ 

HCO+ 1-0 

7+1 

8.35+2 

- 

HNC 1-0 

7+1 

6+2 

- 

HC3N 10-9 

5+1 

6+2 

- 

HCN 3-2 

41+4 

43.3+0.3 

- 

HCO+ 3-2 

24+2 

21+0.3 

- 

HNC 3-2 

39+4 

43+0.3 

37+4 

HC3N 30-29 

15+2 

- 

11 + 1 


(1) This work ; (2) Single-dish fluxes observed with the IRAM 30m telescope ; (3) SMA observations bv \Costaeliola er^il20jJl) . extended conflguration, beam size 
of 0.4". 
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